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PREFACE 


This  report  describes  the  work  completed  during  the  second  year  of  a  two- 
year  study  of  ways  to  enhance  the  use  of  data  from  the  USAF  Defense  Meteor¬ 
ology  Satellite  Program  Topside  Ionospheric  Plasma  Monitor  (SSIE)  at  the  Air 
Force  Global  Weather  Central  (AFGWC).  This  project  is  a  continuation  of  work 
completed  under  contract  to  Boston  College  (Boston  College  Subcontract  No. 
930-1  to  Air  Force  Geophysics  Laboratory  Contract  F19628-82-K-0011)*.  I  wish 
to  acknowledge  the  assistance  of  Dr.  Fred  Rich,  AFGL/PHG;  Ms.  Susan  Bredesen, 
AFGL/PHG;  and  Mr.  8ob  Bussey,  AFGWC/SDDE,  during  this  project  year. 


TABLE  OF  CONTENTS  Page 

1.  INTRODUCTION  j 

2.  PROJECT  DEFINITION  3 

3.  TOPSIDE  PROFILE  MODEL  STUDY  5 

3.1  Task  2:  I R 1 79  Topside  Parametr ization  5 

3.2  Task  6:  Alternative  Topside  Parametrizations  5 

3.2.1  Prof i le  Model  6 

3. 2. 1.1  Model  Development  and  Definition  6 

3. 2. 1.2  Fitting  Profiles  to  Data  11 

3. 2. 1.3  Model  Limitations  20 

3.2.2  Transition  Height  Model  30 

4.  IONOSPHERIC  DATA  PREPROCESSORS  STUDY  37 

4.1  Task  4:  Ne'840)  Data  Preprocessor  Improvements  38 

4.1.1  Spectral  Ne(840)  Model  40 

4.1.2  Effective  F 10  Analysis  41 

4.1.3  Global  Analysis  44 

4.1.4  Preprocessor  Program  (GRIDNE)  Description  45 

4.1.5  Preprocessor  Program  Tests  49 

4. 1.5.1  Description  of  Test  Data  and  Procedures  49 

4. 1.5. 2  Test  Results  54 

4.1.6  Limitations  72 

4.2  Task  5:  Combined  TEC  and  Ne(840)  Preprocessor  73 

r.  4D  MODEL  IMPROVEMENTS  STUDY  80 

5.1  Task  1:  40  Model  Global  Analysis  80 

5.2  Task  3:  Model-based  4D  Height  Functions  80 

6.  PROJECT  SUMMARY  81 

REFERENCES  83 

APPENDIX  A.  Topside  EDP  Parametrizations  85 


iv 


LIST  OF  FIGURES 

Figure  Caption  Page 

1.  SSIE  data  processing  flow  at  AFGWC.  2 

2.  Variation  of  TEC  and  N  (840)  with  various  profile  model  para-  15 

meters.  e 

3.  Results  of  fitting  profiles  to  TEC  and  Ng(840).  19 

4.  The  two  data  sets  used  in  the  TEC  fit  tests.  21 

5.  Variation  of  parameters  used  to  adjust  the  profile  model  to  22 

fit  observed  TEC  for  the  02-06  September  data  set. 

6.  Same  as  Figure  5  for  the  15-19  August  data  set.  24 

7.  Electron  density  profiles  before  (solid  curve)  and  after  26 

(dashed  curves)  adjusting  the  profile  to  fit  TEC  for  the 

1530LT  data  from  the  02-06  September  data  set. 

8.  Same  as  Figure  7  for  the  2030LT  data  from  the  02-06  September  27 

data  set. 

9.  Variation  of  TEC  and  N  (840)  as  functions  of  the  0+/H+  tran-  29 

sition  height  for  fourh/alues  of  yt. 

10.  Seasonal  variation  in  the  h^.  model  for  F10. 7=100,  Kp=3  .  35 

11.  F10.7  variation  in  the  h^  model  for  June  solstice,  Kp=3  .  36 

12.  K  variation  in  the  hT  model  for  June  solstice,  F10. 7=100.  36 

P  T 

13.  Northern  geographic  hemisphere  fields  of  N  (840)  from  the  DHR  12 

(upper  panel)  and  Bent  (lower  panel)  mode1ecoeff  idents  for  a 
sunspot  number  of  120. 

14.  Program  GRIDNE  structure  chart.  46 

15.  Summary  of  the  RMS  analysis  error  between  the  initial,  effec-  52 

tive  F10  analysis,  and  global  analysis  fields  and  the  Bent 

field  for  the  entire  24  hour  analysis  period. 

16.  The  OOOOUT  analysis  grids,  N  (340)  grids,  and  Bent  grids  for  53 

SSN=120,  2  satellites.  * 

17.  The  RMS  difference  between  the  input  f  data  and  f  calculated  55 

from  Equation  22  as  a  function  of  gl8bal  analysfs  iteration 
number. 


v 


56 


IS.  Analysis  results  for  the  first  full  orbit  of  both  DMSP  satel- 
1 i tes  for  test  case  If. 

19.  Summary  of  the  N  (840)  grid  set  comparisons  for  test  case  la  59 

($SN=10,  1  satellite). 

20.  Analysis  results  for  the  0000UT  grid  for  test  case  la  (SSN=10,  60 

1  satellite). 

21.  Summary  of  the  N  (840)  grid  set  comparisons  for  test  case  lb  61 

(SSN«10,  2  satellites). 

22.  Analysis  results  for  the  0000UT  grid  for  test  case  lb  (SSN=10,  62 

2  satel 1 i tes ) . 

23.  Summary  of  the  N  (840)  grid  set  comparisons  for  test  case  lc  63 

(SSN=60,  1  satellite). 

24.  Analysis  results  for  the  0000UT  grid  for  test  case  lc  (SSN=60,  64 

1  satellite). 

25.  Sumnary  of  the  N  (840)  grid  set  comparisons  for  test  case  Id  65 

(SSN=60,  2  satel  fites). 

26.  Analysis  results  for  the  OOOOllT  grid  set  for  test  case  Id  66 

(SSN=60,  2  satellites). 

27.  Su'nnary  of  the  N  (840)  grid  set  comparison  for  test  case  le  67 

( $SN=120,  1  satellite). 

23.  Analysis  restuls  for  the  OOOOUT  grid  for  test  case  le  68 

( SSN=120,  1  satellite). 

29.  Sumnary  of  the  N  (840)  grid  set  comparisons  for  test  case  If  69 

:SSN*120,  2  satellites). 

30.  Analysis  results  for  the  OOOOUT  grid  for  test  case  If  70 

(SSN*120,  2  satellites). 

31.  The  OOOOUT  analysis  and  N  (840)  grids  for  the  Stations,  Grid  76 

1,  and  Grid  2  analyses,  and  the  Bent  OOOOUT  Ng(840)  grid. 

32.  The  OOOOUT  analysis  grids  and  N  (340)  grids  for  the  Stations 

(NT'  and  Grid  3  (NT)  analyses ,  eand  the  Bent  OOOOUT  N  (840) 
grid.  e 


78 


LIST  OF  TABLES 


Caption 

Initial  values  of  the  normalization  constants. 


Page 

IS 


Table 

1. 


2. 

Profile  parameter  values. 

20 

3. 

Constants  for  use  in  Equation 
h.j.  with  solar  epoc^. 

21  describing  the  variation  of 

33 

4. 

Description  of  program  GRIDNE 

routines. 

46 

5. 

Test  parameters. 

5C 

6. 

Results  of  the  effective  F10  analyse'  in  Test  1. 

54 

7. 

Surmary  of  the  results  of  the 

Ng(840)  grid  set  comparisons. 

57 

8. 

Test  2  results  summary. 

/  1 

9. 

Test  3  results  summary. 

~y  o 

10. 

TEC  station  penetration  point 

locatons. 

74 

11. 

TEC  tests  results  summary. 

74 

12. 

Results  of  the  SSIE  effective 

F10  analyses. 

79 

I' <  < 


1.  INTRODUCTION 


Defense  Meteorology  Satellite  Program  (DMSP)  Block  50  satellites  have 

2  3 

Included  Topside  Ionospheric  Plasma  Monitor  (SSIE)  sensors  since  mid  1977  ’  . 
Ionospheric  electron  density,  Ne(840),  and  plasma  scale  height,  8^(840),  from 
the  SSIE  sensors  are  routinely  processed  at  the  Air  Force  Global  Weather 
Center  (AFGWC)  for  specification  of  the  state  of  the  earth's  ionosphere. 
Figure  1  is  a  simplified  flow  diagram  of  data  from  the  SSIE  sensors  within 
AFGWC.  SSIE  data  are  stripped  from  the  DMSP  satellite  telemetry  stream  and 
written  to  a  temporary  file  (IEPREPFILE)  by  the  AFGWC  DMSP  processing  system. 
Program  SSIE  reads  the  raw  SSIE  data  from  this  file,  calculates  electron 

o  3 

densities  and  plasma  scale  heights  ’  ,  and  writes  out  a  data  record  for  every 
64  seconds  containing  an  average  Ng(840)  value  and  Hp(340)  value  to  a 
temporary  *Me  (IFlUXFILE).  Program  LDIPIE  reads  the  data  records  from  this 
file,  reformats  them,  and  writes  them  into  the  Astrogeophysical  Data  Base 
(AGDB)^.  Program  SSIELD  reads  SSIE  Ne(840)  data  for  the  most  recent  24  hours 
from  the  AGDB,  constructs  10°  latitude  by  15°  longitude  northern-hemi sphere 
grids  of  Ng(840)  for  each  of  the  24  hours,  and  writes  them  to  a  gridded-data 
file  (FOURDGRIDS)^.  These  data  are  merged  with  data  grids  from  the  fQF2  and 
TEC  preprocessors  and  written  to  temporary  file  18  by  program  GRDOUT  as 
observation  data  sets2.  The  Air  Weather  Service  (AWS)  4D  Ionospheric  Analysis 
Model**’7  inputs  the  observation  data  sets  from  temporary  file  18  and  produces 
a  consistent  analysis  of  the  ionospheric  electron  density  distribution  over 
the  entire  northern  hemisphere. 

The  primary  objective  of  this  project  is  to  study  methods  for  improving 
the  use  of  the  SSIE  Ng(840)  and  H  (840)  data  by  the  4D  model  system,  identi¬ 
fied  in  Figure  1  by  those  files  and  programs  within  the  dashed  box.  The 
studies  have  focused  on  possible  improvements  to  the  4D  model,  particularly  to 
its  internal  electron-density  profile  model,  and  to  the  various  4D  model 
system  preprocessors. 
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Figure  1.  SSIE  data  processing  flow  at  AFGWC. 

2 


2.  PROJECT  DEFINITION 


The  first  step  in  the  project  was  to  define  specif::  tasks  to  pursue 
within  the  areas  of  research  defined  in  the  project  proposal  and  statement  of 
work.  This  study  resulted  in  the  following  list  of  tasks  to  be  completed 
during  the  project: 

a.  Task  1:  Determine  the  modifications  required  toco”iigure  the  4D  iono¬ 
spheric  analysis  model  for  a  full  global  analysis. 

b.  Task  2:  Investigate  the  use  of  the  profile  parame' r  Ization  used  in  the 

8 

I R 1 79  reference  ionosphere  as  a  replacement  for  the  p3rametrization  cur¬ 
rently  used  in  ‘.he  4D  model. 

c.  Task  3‘  Derive  sets  of  height-profile  basis  func'.ions  for  the  4D  model 
using  plasma  ' requency  profiles  from  the  IRI79  reference  ionosphere. 

d.  Task  4.  Investigate  techniques  for  improving  the  preparation  of  data 
grids  of  Nf;840)  from  the  SSIE  N&(840)  observations. 

e.  T;sk  5:  Work  with  personnel  from  AFGL  and  AFGWC  to  develop  improved 
techniques  for  combining  TEC  and  SSIE  Ne( 840)  data  into  a  single  specification 
of  the  topside  ionosphere. 

f.  Task  6:  Investigate  alternative  parametrizations  of  the  topside  iono¬ 
sphere  which  would  use  the  full  set  of  data  available  from  the  SSIE  sensors 
and  all  other  available  data,  including  Total  Electron  Content  (TEC)  data,  to 
specify  the  electron  density  profile  above  the  F2-layer  peak. 

g.  Task  7:  8egin  an  investigation  into  ways  to  use  data  from  the  SSIE 
sensors,  in  conjunction  witn  other  observations,  for  modeling  the  ionospheric 
subauroral  trough  as  part  of  the  4D  model  preprocessor  system. 

These  tasks  fall  generally  into  three  main  studies:  1)  Topside  Profile 

Model  Study  (Tasks  2  and  6),  2)  Ionospheric  Data  Preprocessor  Improvements 

Study  (Tasks  4,  5,  and  7),  and  3)  4D  Model  Improvements  Study  (Tasks  1  and  3). 

Tasks  1,  2,  and  3  were  completed  during  the  first  year  of  the  study  and  the 

g 

outcomes  reported  in  Scientific  Report  No.  1  of  this  contract  .  Task  7  was  to 
be  worked  as  resources  allowed.  As  the  work  on  improvements  to  the  N^JSAO) 


gridding  program  (Task  4)  grew  in  scope  during  the  course  of  the  project, 
there  were  no  resources  available  to  put  on  Task  7.  The  work  completed  in 
Tasks  4,  5,  and  6  will  be  reported  here.. 


3.  TOPSIDE  PROFILE  MODEL  STUDY 


This  study  includes  Tasks  2  and  4  of  the  project.  The  overall  objective 
of  this  study  4s  to  investigate  improved  parametri zations  of  the  topside 
Ionospheric  electron  density  profile  for  use  within  the  40  model. 

3.1  Task  2:  I R 1 79  Topside  Parametrization 

The  major  thrust  of  this  task  was  to  investigate  the  I R 1 79  model  to 

determine  if  1)  the  I R 1 79  topside  parametrization  could  be  modified  for  use 

within  the  4D  model,  and  2)  the  profile  shape  provided  by  the  I R 1 79  was 

sufficiently  better  than  the  Damon-Hartranft-kamsay  (DHR)  profile  model^’** 

currently  used  by  the  40  model  to  justify  replacing  the  DHR  model  with  I R 1 79 . 

q 

As  reported  earlier  ,  the  conclusion  of  the  investigation  was  that  the  I R 1 79 
was  not  a  viable  replacement  for  the  DHR  model,  so  work  was  terminated  on  this 
task  and  begun  on  Task  6. 

3.2  Task  6:  Alternative  Topside  Paramatrizations 

The  purpose  of  this  task  was  to  investigate  parametr izations  of  the 
topside  electron  density  profile  (EDP)  wh’ch  would  represent  a  two-consti¬ 
tuent  (0+  and  H+ )  topside  ionosphere  and  which  could  be  easily  adjusted  to  fit 
observed  conditions.  The  objective  was  to  find  a  parametrization  which  could 
use  the  full  complement  of  observations  available  from  the  SSIE  (and  SSIESl 
sensor,  which  includes  the  0+  and  H+  densities,  N(0  )  and  N(H  ),  and  the  elec¬ 
tron  and  ion  temperatures ,  T^  and  T^,  in  conjunction  with  ground-based  obser¬ 
vations  of  f  F2,  h  F2,  and  TEC.  It  was  decided  during  the  first  year  of  the 
o  ’  m  +  + 

project  to  base  the  parametrization  on  a  model  of  the  0  to  H  transition 
height,  hy,  defined  as  the  height  at  which  N(0+)  3  N(H+).  This  task  then 
divided  into  two  subtasks:  an  investigation  of  EDP  models  based  on  hT,  and  the 
development  of  a  model  for  h^.. 
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3.2.1  Profile  Model 


3. 2. 1.1  Model  Development  and  Definition 

Several  different  topside  electron  density  profile  (EDP)  model  parame- 
trizations  were  investigated  as  part  of  this  task.  These  included  modified 
versions  of  the  AWS  DHR  model  and  the  Bent  model  ,  two  models  based  on  a 
parametrization  of  the  topside  scale  height  profile  developed  to  analyze 
topside  sounder  data^,  and  a  parametrization  based  on  a  two-component  diffu¬ 
sive  equilibrium,  (DE)  model.  For  the  purpose  of  calculating  TEC,  all  models 
used  the  DHR  bottomside  model  for  the  EDP  below  the  F2  peak.  These  models 
were  evaluated  against  the  following  criteria: 

1.  The  electron  density  variation  with  height,  N  (h),  should  be  con¬ 
tinuous  and  Smooth. 

2.  The  scale  height  variation  with  height,  H(h),  should  be  continuous 
and  reasonably  smooth  (i.e.,  as  few  sharp  discontinuities  as  possible). 

3.  The  model  should  be  capable  of  fitting  a  wide  range  of  data  types, 
including  fQF2,  hmF2,  TEC,  and  the  SSIE  data  set,  but  should  require  only  hmF2 

and  f  F2. 

0 

4.  Conversely,  the  model  should  reflect  the  basic  physical  processes 
active  in  the  topside  ionosphere  such  that  it  will  not  produce  physically 
unreasonable  profiles  in  attempting  to  fit  an  inconsistent  data  set.  This 
would  be  useful  in  data  quality  control  and  in  intercomparison  of  data  sets. 

5.  As  implemented  on  a  computer,  the  profile  model  and  the  algorithms 
used  to  fit  the  profile  to  observations  should  be  as  simple  and  computa¬ 
tionally  quick  as  possible  and  computationally  stable. 

The  first  four  parametri zat ions  considered,  which  are  briefly  described 
in  Appendix  B,  all  met  the  first  two  criteria,  but  ran  into  difficulties  with 
one  or  more  of  the  last  three,  particularly  numbers  three  and  four.  The 
primary  reasons  for  rejecting  the  various  models  were  as  follows: 

1.  Modified  DHR  model.  Although  initially  attrocti/e,  this  model  was 
abandoned  fairly  early  for  a  number  of  reasons.  In  order  to  produce  a  profile 


with  only  hmF2  and  fQF2  input,  a  model  would  be  required  to  provide  the  scale 
height  for  the  exponential  layer  above  h=hy,  and  the  DHR  scale  height  model ^ 
would  need  to  be  modified  to  account  for  the  addition  of  the  exponential 
layer.  Both  actions  would  require  analysis  of  large  amounts  of  data  which 
were  not  available.  Additionally,  the  procedures  to  fit  the  profile  to 
observed  SSIE  data  became  very  cumbersome,  and  could  very  easily  produce 
unreal istic  profiles. 

2.  Modified  8ent  model.  This  parametruation  was  developed  as  a  solu¬ 
tion  to  the  1  imi ted-data-prof i le  problem  encountered  with  the  modified  DHR 
model.  Unfortunately,  the  data-fitting  procedures  proved  to  be  as  cumbersome 
as  with  the  modified  DHR  model,  and,  in  some  cases,  the  model  produced  N  ( h ) 
curves  which  increased  with  height. 

3.  Titheridge  H(h)  model.  This  model  was  considered  in  an  attempt  to 
control  the  physically  unrealistic  behavior  exhibited  by  the  first  two 
models.  As  with  the  first  two  models,  the  data-fitting  procedures  were  a 
problem.  Not  only  were  they  cumbersome,  but  they  were  also  computationally 
unstable  under  reasonable  conditions.  Although  the  Titheridge  model  was 
difficult  to  use  in  fitting  observations,  this  was  more  a  problem  with  the 
particular  parametrization  chosen  than  an  inherent  problem  with  a  diffusive 
equilibrium  (DE)  distribution.  This  led  to  developing  a  parametrization  of  a 
DE  profile  more  suited  to  the  present  application. 


The  height  distribution  of  the  ion  density  for  a  topside  ionosphere  in 
diffusive  equilibrium  is  given  by  ^ 
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where  N  (x.)  is  the  density  of  ion  species  j  at  the  reference  height  h  ,  T  -is  the 

0  J  o  o 

plasma  temperature  (Tg  +  T.),  ok  is  the  mass  of  species  j,  and  m  is  the  mean 
ionic  mass  given  by 

m  =  ( 7:  N(X 

'  J 

Collecting  the  terms  which  are  species-independent,  Equation  1  can  be 
rewritten 


where 


dh 


h 


Integrating  the  second  term  in  Equation  2  then  gives 

h 

Ni  s  e  £  VXj}  exp  (‘f  kT7  dh) 

P  j  v  1  ' 

where 


(2) 


For  an  ionosphere  consisting  of  0+  and  H  ions  in  diffusive  equilibrium,  the 
mean  ionic  mass  can  be  simplified  to 


m  =  m 


16+R 

UR 


where  m+  is  the  mass  of  H+,  and  R  is  the  ratio  of  N(H+)  to  N(0+)  given  by 
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r  s  NiiLi  =  r  exp 
N(0+)  0 


/ 

hi-i 


dh 


(3) 


Thus,  for  an  0+,  H+  ionosphere  in  diffusive  equilibrium,  the  density  distri¬ 


bution  becomes 


T 

rf 

e 

No(0+)  e’161  ♦  Nq(H+)  e_I 

l  r>  J 

(4) 


dh 


where 


(5) 


This  provides  a  topside  model  which  requires  a  model  for  T  (h)  and  T.(h),  some 

method  for  determining  Nq(0+)  and  Nq(H+),  and  a  way  to  couple  the  profile  to 

the  F2  peak  at  h  F2. 

m 

A  recent  empirical  model  based  on  data  from  the  AE-C  satellite^  was 
chosen  to  describe  the  height  variation  of  T^  and  T^.  In  this  model,  both 
parameters  are  represented  by  equations  of  the  form 


where  the  coefficients,  a  are  functions  of  geomagnetic  latitude  and  local 

time.  In  the  original  model,  summer  solstice  values  for  T . ( h )  were  given  for 

J 

northern  magnetic  latitudes  and  winter  solstice  values  were  given  for 
southern  magnetic  latitudes.  This  was  modified  to  provide  a  seasonal  varia¬ 
tion  at  all  latitudes  by  multiplying  the  non-symmetr ic  (in  latitude)  terms  in 
the  expression  describing  the  anj  coefficients  by  cos  [ ( 2  t/ 3 6 5 )  (D  +  192)] 
where  D  is  the  day  of  the  year. 


It  was  decided  to  use  the  Bent  topside  parabolic  layer  to  describe  the 
profile  from  hmF2  to  the  base  of  the  OE  model  for  three  reasons.  First,  a 
parabolic  layer  approximates  the  shape  of  the  EDP  near  the  F2  peak  fairly 
well;  second,  it  is  a  mathematically  simple  representation,  with  the  density 
and  scale  height  given  by 


(7) 
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A  -  (h-v?>2 

h-h  F? 
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where  yt  is  the  semi-thickness;  and  third,  the  Bent  model  values  for  y^  can  be 
used  in  conjunction  with  a  model  for  the  0+  to  H+  transition  height,  hy,  to 
specify  Nq(0+)  and  Nq(H+),  thereby  requiring  only  fQ F2  and  h^F2  to  specify  the 
entire  profile.  This  is  done  by  determining  the  height,  h^,  at  which  the 
scale  height  given  by  Equation  8  is  equal  to  that  calculated  from 


-1 


(9) 


where  N  (h)  is  given  by  Equation  4  (assuming  N  (h)  *  N.(h)),  and  using  the 
e  ®  Y  + 

density  at  that  height  from  Equation  7  to  determine  Nq(0  )  and  Nq(H  ). 

It  is  possible  to  do  this  because  the  scale  height  defined  by  Equations 
9  and  4  is  not  a  function  of  N(0+)  or  N(H+),  but  rather  of  their  ratio,  R. 
For  an  0+/H+  ionosphere,  Equation  1  can  be  rewritten 


Ni  =  N(0+)  +  N(H+) 


(10) 


where 


h 


Combing  Equations  9  and  10, 

H  , _ _ 

N(0+)f(0+)  ♦  N(H+)f(H+) 

or,  from  Equation  3, 
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DE  ' 


_ Xjl  A _ 

f(0+)  ♦  Rf(H+) 


(ID 


10 


Using  the  expression  for  f(X+)  from  Equation  11,  this  becomes 

<3 


(121 


Since  m  is  only  a  function  of  R,  is  not  a  function  of  the  densities,  Nq( 0+ ' 

and  N  (H+),  but  rather  their  ratio,  R„. 
o  o 

Given  values  for  h  F2,  f  F2,  and  hT,  this  fact  can  be  used  to  fit  the 

m  o  T 

parabolic  and  DE  profiles  together  as  follows: 

1.  Set  hQ  to  h-j..  Since  hv  is  defined  as  the  height  where  N(0+)  = 
N(H+),  then  Rq  -  1.0. 

2.  Using  Equations  8  and  12  with  the  temperature  model  (Equation  6', 
iteratively  locate  the  height,  hj,  at  which  the  scale  heights  from  Equations  S 
and  12  are  the  same. 

3.  Calculate  the  total  ion  (electron)  density  at  this  height,  N  (h.), 

+  +  G  1 

from  Equation  7.  The  densities  for  0  and  H  at  h^  can  then  be  calculated  from 
+  N  (h  ) 

Ni(0  }  s  •  and 

N1(H+)  =  R(h1)N1(0+). 

4.  Set  hQ  to  hj,  RQ  to  R(h1),  Nq(0+)  to  Nj(0+).  and  Nq(H+)  to  N^H*). 
All  profile  parameters  are  now  calculated. 

3. 2. 1.2  Fitting  Profiles  to  Oata 

The  various  types  of  inospheric  observations  available  are  used  to 
adjust  the  profile  model  as  follows: 

a.  The  ion  densities,  N^(o+)  and  NS(H+),  the  ion  temperature,  T^,  and 
the  electron  temperature,  Tes,  at  the  satellite  altitude  are  used  to  calculate 
the  0  /H  transition  height,  h^. 
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b.  The  temperatures,  Tis  and  Tg$,  can  be  used  to  adjust  the  AE-C 
temperature  model . 

c.  The  electron  density  at  the  satellite  (Ne(840)),  fQF2,  hmF2,  and  ^ 
are  used  to  adjust  various  model  parameters  so  that  the  resulting  profile 
matches  the  observations. 

As  stated  earlier,  the  only  data  required  in  order  to  build  a  profile  are  fQF2 

and  h  1-2.  All  other  observations  are  used  as  available, 
m 

The  transition  height  is  calculated  from  Equation  3  for  the  height  varia¬ 
tion  of  the  ratio  of  N(H+)  to  N(0+).  Assuming  that  g  is  constant,  this 
equation  becomes 


at  the  transition  height,  and 


at  the  satellite  altitude,  where  G  r  15in+g/k.  Combing  these  yields 


840 


If  T  is  assumed  to  vary  linearly  with  height  between  hy  and  840km,  the 
integral  becomes 


84  j 


dh _ 
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so  that 
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c 
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S 

Ml 


T. 

h  =  840  +  4i  (Rj  -  1)  (13) 

-  i  s 

where  i  -  ^T.k/15m+g.  In  order  to  reduce  errors  introduced  by  assuming 
constant  values  for  g  and  the  transition  height  is  calculated  iteratively 
as  follows: 

1.  Set  g  to  its  value  at  840km,  calculate  .IT.  at  840km  from  the  AE-C 
model,  and  calculate  an  initial  estimate  for  hy  from  Equation  13. 

2.  Calculate  a  new  estimate  for  hy  using  values  for  g  and  :.T.  at  a 
height  midway  between  the  last  estimate  for  hy  and  840km. 

3.  Iterate  on  step  2  until  the  change  in  hy  from  one  iteration  to  the 
next  is  below  some  threshold. 

Values  for  hy  calculated  in  this  fashion  can  be  used  either  in  lieu  of  the  hy 
model  described  in  section  3.2.2  or  to  update  and  modify  this  model. 

Modifying  the  profile  mcdel  to  fit  observed  <qF2  and  h^?  with  either 
Ne( 840)  or  TEC,  or  both  Ng(840)  and  TEC,  is  also  an  iterative  process.  Three 
profile  adjustment  parameters  were  added  to  the  equations  describing  the 
profile  to  facilitate  the  fitting  process.  The  basic  change  is  to  Equation  4, 
which  becomes 


Ni  * 
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No(0+)  e'16*1  ♦  N0(H+)e*1 


‘0  +  ^  (h-400), 

where  8 ,  a  ,  and  ^  are  the  adjustment  parameters.  This  change  requires 
modification  of  Equations  5  and  12  to 


R  =  R^16'1'1*1 


and 
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where 
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This  provides  a  total  of  four  model  parameters  which  can  be  adjusted  In  order 
to  fit  a  profile  to  observations  -  a  ,  S,  and  y^. 

The  variation  of  Ng(840)  and  TEC  as  functions  of  these  parameters  is 

illustrated  in  Figure  2.  All  calculations  were  for  mid-latitude  at  local  noon 

with  f  F2  =  7  MHz  and  h  F2  -  350km,  and  the  four  curves  in  each  plot  are  for 
o  m 

fojr  values  of  h-j.:  650kmn,  1000km,  1500km,  2000km.  In  all  plots,  the  top 
curve  is  for  hy  =  650km  and  the  bottom  is  for  hy  *  2000km.  Also  included  in 
this  figure  are  plots  showing  the  variation  of  Ng(840)  and  TEC  as  functions  of 
f„F2  and  T  (400),  the  plasma  temperature  at  400km,  for  comparison. 

The  procedure  developed  for  adjusting  the  profile  parameters  to  fit 
observed  TEC  and  Ne(840)  allows  for  a  controlled  iterative  adjustment  to  any 
three  profile  parameters  until  an  adequate  fit  is  achieved.  This  procedure 
involves  the  solution  to  the  equations 
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(14) 

(15) 


under  the  contraint  that  a  normalized  change  distance,  Al  ,  given  by 


-L 


(16) 


is  minimized.  In  these  equations,  T  and  N  designate  TEC  and  Ng(840);  (x^, 
i=l,3)  are  any  three  profile  parameters;  and  (c^»  1=1,3)  are  three  normaliza¬ 
tion  factors.  This  system  of  equations  is  solved  by  reducing  Equation  16  to  a 
single  tx.  variable  using  Equations  14  and  15,  differentiating  this  with 
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respect  to  ax^,  and  solving  for  Ax^.  This  yields 
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If  only  one  of  TEC  or  Ne(840)  is  available,  a  similar  set  of  equations  can  be 
used  to  adjust  any  two  parameters  for  the  fit.  For  TEC,  these  are 

A  .i  Ti  T 
Ax, - X  -n  at  , 

d  4  (18) 
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where 
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Initially,  the  normal ization  constants  were  included  to  account  for  the 
different  ranges  of  the  various  parameters  and  were  set  to  the  values  shown  in 
Table  1.  It  was  found,  however,  that  the  amount,  x^,  by  which  a  parameter 
would  be  varied  could  be  controlled  by  adjusting  its  c^  value  with  respect  to 
the  c.  values  of  the  other  parameters  being  modified. 


Table  1.  Initial  values  of  the  normalization  constants. 


Parameter 


c.  value 
i 


al 

0.0004 

0.07 

~0 

0.07 

20.0 

V2 

1.0 

Tp(400) 

200.0 

hmF2 

20.0 

hT 

100.0 

As  there  were  no  reliable  sets  of  Ne(840)  observations  available,  only  a 

few  limited  tests  were  run  of  the  procedures  developed  for  fitting  profiles  to 

both  TEC  and  Ne(840)  data.  Figure  3  illustrates  the  results  of  one  of  these 

tests.  The  solid  curve  is  a  profile  calculated  for  an  f  F2  of  8.4  and  h  F2  of 

am 

314km  for  a  location  of  20°N,  203°W  (Palehua  Observatory  in  Havaii)  at  1930 

16  2 

local  time.  This  profile  has  a  TEC  of  21.8  x  10  el/m  and  an  N  (840)  of  4.28 

4  3  e 

x  10  el/cm  .  The  fitting  procedure  using  Equations  17  was  then  used  to  fit  to 

the  same  Ng(840)  for  different  values  of  TEC  by  adjusting  parameters  y^, 
and  k  .  The  upper  and  lower  dashed  curves  in  Figure  3  are  the  profiles 
produced  for  TEC  values  of  28.0  x  10*62l/in2  and  20.0  x  10^el/m2,  respec¬ 
tively.  Table  2  lists  the  values  of  the  three  adjusted  parameters  for  each 
profile.  All  tests  showed  this  to  be  a  very  flexible  and  robust  technique. 
Further  testing  and  evaluation  of  this  technique  must  await  the  availability 
of  reliable  in  situ  data  from  the  SSIES  sensor. 
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Figure  3.  Results  of  fitting  profiles  to  TFC  end  N  '840). 


Table  2.  Profile  parameter  values. 


II c 

h 

p 

ao 

20.0 

69.0 

0.985 

0.938 

21.8 

136.0 

1.000 

1.000 

28.0 

296.0 

0.942 

1.318 

(xl016el/m2) 

(km) 

A  more  extensive  series  of  tests  was  carried  out  on  the  procedure  for 

fitting  profiles  to  f  F2,  h  F2,  and  TEC  data.  A  set  of  near-coincident 

o  m 

observations  of  TEC  from  the  AWS  Palehua  Observatory  and  observations  of  fQF2 
and  h  F2  (calculated  from  the  M3000  factor)  from  an  lonosonde  on  Maul  taken 
during  August  and  September  1983  were  used  to  Investigate  the  fitting  proce¬ 
dure  using  Equations  18.  Two  twenty-four  hour,  five-day  mean,  data  sets  were 
constructed  for  two  geomagnetical 1y  quiet  periods  -  one  centered  on  17  August, 
the  other  centered  on  4  September.  Figure  4  shows  the  f  F2,  h  F2,  and  TEC 
variations  for  these  two  data  sets.  The  dashed  curve  In  the  two  TEC  plots  Is 
the  TEC  calculated  from  the  OE  profile  model  using  only  the  fQF2  and  hmF2 
data.  Figures  5  and  6  show  the  adjustments  required  to  fit  the  observed  TEC 
(within  S*)  for  four  pairs  of  adjustment  parameters  -  and  clq,  6  and  aQ,  y ^ 
and  .  ,  and  ;  and  y, .  Also  shown  are  two  cases  In  which  f  F2  was  allowed  to 
vary  only  by  small  increments  by  settng  c^  for  f  F 2  to  0.3.  In  all  cases*  the 
TEC  fit  was  achieved  in  less  than  5  interations,  with  2-3  iterations  being 
typical.  To  illustrate  the  effects  of  the  various  modifications  on  the  entire 
profile.  Figures  7  and  8  show  the  profiles  for  1530  IT  and  2030  IT  from  the  02- 
05  September  data  set  which  represent  the  maximum  negative  and  positive  TEC 
corrections  required,  respectively.  The  solid  curve  shows  the  Initial  pro¬ 
file,  and  the  dashed  curves  show  profiles  which  match  the  observed  TEC  through 
adjustment  of  the  indicated  parameters. 

3.2. 1.3  Model  Limitations 

Many  of  the  potential  limitations  of  this  profile  model  derive  from  the 
assumptions  made  in  developing  Equation  4.  The  most  basic  of  these,  that  the 
entire  profile  is  in  diffusive  equilibrium,  is  not  too  bad  an  assumption  a 
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good  pa^-t  of  the  time.  In  a  study  of  departures  from  diffusive  equilibrium  in 
the  topside  ionosphere,  Lockwood  and  Titheridge^  found  that  roughly  37X  of  the 
more  than  10,000  topside  soundings  they  analyzed  showed  appreciable  depar¬ 
tures  from  diffusive  equilibrium.  The  most  serious  of  these  were  in  the 
equatorial  regions  (c.f.  their  Figure  8)  where  large  0+  fluxes  can  occur  which 
can  seriously  effect  the  shape  of  the  profile  just  above  the  F2  peak  and  can 
cause  appreciable  errors  in  TEC  calculated  by  the  model.  Additionally,  large 
vertical  fluxes  of  0+  and  H+  can  be  found  poleward  of  the  plasmapause  location 
wrnch  can  also  seriously  deform  the  vertical  density  profile.  Although 
neither  effect  is  included  in  the  current  model,  a  systenatic  study  of  the 
adjustments  required  to  fit  profiles  to  data  could  lead  to  empirical  models 
for  the  adjustment  parameters,  5,  -v  ,  and  which  would  improve  the  model's 
ability  to  produce  realistic  profiles  with  only  fQF2  and  hmF2 . 

Another  possible  problem  area  is  the  assumption  that  only  0+  and  H+  are 
present  in  the  ionosphere.  While  this  may  have  only  a  small  effect  on  TEC  or 
Ne ( 340 )  calculated  from  the  model,  it  could  have  a  serious  effect  on  the 
procedure  for  calculating  hy  from  the  SSIES  data  set.  The  procedure  assumes 
that  the  light  ion  density  calculated  from  the  ion  RPA  analysis  is  H+,  an 
assumption  that  may  not  always  be  valid.  If  the  density  were  actually  that 
for  He+,  the  estimate  for  hy  would  be  too  low.  Figure  9  illustrates  the 
variation  of  TEC  and  Ne(840)  as  functions  of  hy  for  the  same  conditions  as 
those  used  in  generating  Figure  5  for  various  values  of  y^.  As  can  be  seen, 
the  effect  of  an  error  in  hy  on  TEC  and  Ne(840)  is  very  small  unless  hy  is 
below  about  1000km.  If  hy  is  below  1000km,  which  occurs  curing  tte  night  at 
mi d- 1 ati tudes  and  during  the  day  at  mid-latitudes  during  solar  minimum  winter 
solstice  conditions  (see  Figure  9),  an  error  of  *  100km  in  hy  will  result  in 
roughly  a  +7%  error  in  TEC  and  +12%  error  in  Ng(840).  Note,  however,  that  if 
the  data  are  available  to  calculate  hy,  then  a  value  for  Ne(840)  is  also 
available.  By  adjusting  the  profile  to  fit  Ng(840),  the  error  in  TEC  intro¬ 
duced  by  an  error  in  hy  should  be  reduced  to  a  few  percent  or  less. 

A  th-""d  assumption  made  in  developing  the  OE  model  is  that  the  magnetic 
field  lines  are  nearly  vertical,  an  assumption  which  is  reasonably  valid  to 
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Figure  9.  Variation  of  TFC  and  N  (349)  as  functions  of  the 
transition  height  for  rou r  values  o*  .  . 


within  roughly  t30°  of  the  dip  equator.  This  is  implied  in  writing  Equation  1 
as  integral  in  height  rather  than  in  distance  along  the  magnetic  field.  In 
the  equatorial  region,  the  breakdown  in  this  assumption  will  have  two  effects 
on  the  profile  generated  by  this  model.  First,  the  profile  shape  produced  by 
the  model  will  not  include  distortions  caused  by  the  non-vertical  field  lines, 
an  effect  which  may  be  somewhat  alleviated  by  the  use  of  an  empirical  model 
for  the  0+/H+  transition  height.  Second,  the  fQF2  and  hmF2  used  in  construct¬ 
ing  the  profile  shou^  not  necessarily  be  values  vertically  below  a  given 
point  in  the  profile,  but  rather  values  more  representative  of  conditions  in 
the  F2  region  at  the  base  of  the  field  line  passing  through  the  point.  While 
not  conceptually  a  problem,  constructing  large  numbers  of  vertical  profiles 
based  on  latitudinal  variations  in  the  F2  region  could  be  computationally 
burdensome. 

Finally,  although  the  data-fitting  techniques  described  in  Section 

2.2. 1.2  appeared  to  be  fairly  robust  during  the  limited  testing  conducted,  a 
full  series  of  tests  to  determine  which  parameters  to  adjust  under  which 
conditions  should  be  run  when  reliable  data  become  available  from  the  SSIES 

sensor . 

3.2.2  Transition  Height  Model 

A  global  model  of  the  0+  to  H+  transition  height  was  developed  from 
published  analyses  cf  topside  sounder  data  and  in  situ  retarding  potential 
analyzer  (RPA)  data.  Titheridge^  analyzed  topside  sounding  data  from  the 
Alouette  I  satellite  to  determine  the  transition  height  variation  during 
solar  minimum  (1964-1965)  and  a  limited  amount  of  data  from  the  ISIS  1  and 
ISIS  2  satellites  near  solar  maximum  (1969  and  1971).  Kutiev  et  ai^’^ 
analyzed  RPA  data  from  the  0G0-6  satellite  to  determine  the  (nighttime)  varia¬ 
tion  of  hT  during  solar  max’mum  (1969-1970).  (It  should  be  noted  that  Kutiev 
et  all°  nave  proposed  the  use  of  the  transition  height  as  an  improvement  to  the 
ISI79  model.)  The  latitudinal,  seasonal,  local  time,  and  solar  epoch  varia¬ 
tions  in  the  hj  model  were  derived  from  these  references.  The  location  of  the 

plasmapause  and  its  variation  with  local  time  and  magnetic  activity  were 
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derived  from  a  further  analysis  of  the  Alouette  hT  data  by  Tithe-Mdge  and 
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work  by  Spiro  reported  in  Moffett  and  Quegan  . 


The  functional  form  of  the  model  Is  as  follows: 
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where  hg,  h^,  and  hp  are  the  equatorial  ('=0°),  mid-latitude,  and  polar  values 
of  hy  which  are  functions  of  local  time,  season,  and  solar  epoch;  •  is  the 
magnetic  (F-layer  apex)  latitude;  and  W  and  W  are  feature  widths.  The 

C  P 

second  and  third  terms  In  the  left-hand  side  of  Equation  19  model  a  nighttime 
equatorial  peak  in  the  latitudinal  variation  of  hy  and  the  transition  into  the 
polar  cap  ionosphere,  respecti vely. 


The  local  time  (LT)  variation  of  Iv,  h.,  and  h  is  modeled  usinq  the 

6  *  fll*  P  5 

epstein  step  function  used  to  model  local  time  variations  in  the  I R 1 7 9  model. 
The  form  of  this  variation  is 


f (LT)  *  f 
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where  fN  and  fg  are  the  night  and  day  values  of  the  function,  and  LT<-R  and  tT<.<. 
are  the  local  times  of  sunrise  and  sunset. 


The  seasonal  variation  of  hgj  h^,  and  hp  is  modeled  by  a  simple  cosine 
variation  of  the  form 


f (day)  =  (fw  ♦  fs) 


(f„  -  fs)  COS 


2  n 

36$ 


( D-10) 


where  f  and  f  are  the  winter  and  summer  values  of  the  function,  and  D  is  the 
w  s 

day  of  the  year.  The  plus  (minus)  sign  in  this  equation  is  used  for  the 
northern  (southern)  hemisphere. 

The  solar  epoch  variation  of  hg>  h^,  and  hp  is  modeled  as  a  linear 
variation  with  the  10.7  cm  solar  radio  flux  (F10.7)  of  the  form 
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where  S  =  75  if  F10.  7  75,  S  =  F10.7  if  75  •  F10.7  180,  and  5  *  180  if 

F 10- 7  ^  180.  Taole  3  shows  the  values  of  h.  and  a,  for  the  three  latitudinal 

10  i 

regimes  (equatorial,  mid,  and  polar)  for  day  and  night,  summer  and  winter. 
There  are  no  day  values  for  the  equatorial  regime,  as  the  equatorial  peak  is 
present  only  at  night. 

The  nighttime  equatorial  peak  term  in  Equation  19  is  included  to  model  a 
peak  present  in  the  Alouette  data  {c.f.  Figure  7  in  Titheridge^)  and  to  fit 
the  variations  presented  in  Kutiev  et  al^.  it  is  modeled  as  a  gaussian  in 
latitude  centered  on  the  magnetic  (dip)  equator  with  a  half-width  of  10°. 

The  transition  through  the  plasmapause,  the  third  term  in  Equation  19,  is 
modeled  through  variations  in  two  latitude  boundaries:  the  equatorward 
boundary  of  the  plasmapause  transition  region,  \j,  and  the  center  of  the 
light-ion  trough,  ^ •  The  transition  is  modeled  as  an  error  function  centered 
on  9  (  '3"  ll/2  with  a  transition  width,  The  ■  boundary  varies  only 

with  local  time  using  Equation  20  with  fg  =  30\  f^  =  50°,  LT<.R  =  9,  LT^  *  20. 
This  variation  was  chosen  to  fit  the  variation  shown  in  Figure  5  of 
Titne’-idge  .  The  boundary  varies  with  local  time  and  magnetic  activity 

V  as 
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This  equation  reflects  variations  with  K  as  reported  by  Titheridge  ,  and 

P 

local  time  variation  as  described  by  Moffett  and  Queganc  and  shown  in  their 
Figure  10. 

A  listing  of  tne  Fortran  function  that  implements  this  model  is  given  in 
Appendix  3. 
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Table  3.  Constants  for  use  In  Equation  21  describing  the  variation  of  with 
solar  epoch. 


Reqime 

Summer 

Day 

Winter 

Sunnier 

Night 

Winter 

Polar 

1200 

1200 

1100 

900 

Mid 

900 

600 

600 

450 

Equatorial 

- 

800 

800 

Day  Night 


Reqime 

Summer 

Winter 

Summer 

Winter 

Polar 

450 

450 

550 

750 

Mid 

600 

700 

300 

200 

Equatorial 

• 

• 

300 

300 
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Figures  10-12  illustrate  the  seasonal,  solar  epoch,  and  magnetic  varia¬ 
tions  in  hj  produced  by  this  model.  In  Figure  10,  the  latitudinal  variation 
of  by  at  noon  and  midnight  are  shown  for  the  June  and  December  solstices  and 
for  an  equinox  for  F10.7  =  100  janskys  and  Kp  =  In  Figure  11,  the 
variation  with  F10.7  is  shown  for  noon  and  midnight  at  the  June  solstice,  Kp  » 
3°.  The  variation  with  Kp  is  shown  in  Figure  12  for  noon  and  midnight  at  the 
June  solstice,  F10.7  =  100  janskys. 

The  primary  weaknesses  in  this  model  are  in  the  solar  epoch  variability 
and  at  high  latitudes.  Only  a  small  amount  of  data  was  available  at  solar 
maximum,  and  most  of  it  was  nighttime  data.  It  may  be  possible  to  refine  this 
variation  when  data  from  the  SSIES  RPA  sensor  becomes  available.  At  high 
latitudes,  the  model  provides,  at  best,  only  a  rough  estimate  of  the  location 
of  the  sharp  transition  from  mid-latitude  to  polar  transition  heights,  and 
detailed  structures  such  as  the  effect  of  the  dayside  cusp  (c.f.  Figure  7  In 
Titheridge^)  are  not  included  In  the  model  due  to  a  lack  of  data.  Again, 
this  section  of  the  model  could  be  refined  using  SSIES  RPA  data. 
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Fiyure  11.  F10.7  variation  In  the  h^  model  for  June  solstice,  K  «3o.  The  solid 

curves  are  for  F10. 7*125,  the  upper  dashed  curves  are  for  F 10 . 7*175, 
and  the  lower  dashed  curves  are  for  F 10 . 7*75 . 
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Figure  12.  Kj  variation  in  the  h_  model  for  June  solstice,  F10.7-100.  The  solid 

curves  are  for  k  =3+,  the  upper  dashed  curves  ere  for  K  =>,  and  the 

lower  dashed  curves  are  for  K  aOr>  u 
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4.  IONOSPHERIC  DATA  PREPROCESSORS  STUDY 

Limited  Investigations  Into  analysis  stability  problems  with  the  40 

model  led  to  the  development  of  front-end  data  preprocessors  for  f  F2  and  TEC 

21  0 
data  by  AFGWC/TSIS  .  These  preprocessors  use  simple  curve-fitting  tech¬ 
niques  to  generate  equally  spaced  grids  of  f  F2  and  TEC  values  from  data 

o  22 

available  In  the  Astrogeophyslcal  Data  Base  ( AGOB )  .  In  the  contract  prior 

21 

to  the  current  one,  a  40  model  preprocessor  system  was  designed  around  these 
two  preprocessors,  programs  FOR  I ER  (f  F2  data)  and  POL Yl  (TEC  data),  and  a 
third  preprocessor  for  SSIE  Ng(840)  data,  program  SSIELD  ,  provided  by 
another  contractor. 

This  system,  tested  and  installed  at  AFGWC,  was  to  be  the  starting  point 
for  the  tasks  under  this  study  area.  Unfortunately,  as  of  the  start  of  this 
phase  of  the  project,  the  SSIELD  program  was  still  not  operational  at  AFGWC, 
and  was  neither  sufficiently  tested  nor  adequately  documented  to  be  of  any 
use.  Thus,  this  study  area  was  begun  by  assisting  AFGL/PHG  with  an  assessment 
of  the  SSIELD  program. 

g 

The  assessment,  documented  In  Scientific  Report  No.  1  ,  showed  that  the 
SSIELD  program  was  not  working  correctly  and,  as  currently  designed,  would  be 
Inadequate  for  the  task  even  If  the  problems  In  the  program  were  corrected. 
The  assessment  results  were  discussed  with  Dr.  Frederick  Rich  (AFGL/PHG)  and 
Mr.  Bob  Bussey  (AWS/SDDE),  and  It  was  decided  to  redirect  the  efforts  in  this 
study  area  to  identifying  and  Investigating  alternative  analysis  techniques 
for  producing  gridded  fields  of  SSIE  Ng(840)  data. 

Toward  this  end,  Informal  discussions  were  held  with  Dr.  Rich  and  Mr. 
Bussey  to  determine  the  requirements  to  be  met  by  the  SSIELD  program.  A 
distillation  of  these  discussions  provided  the  following  list  of  requirements 
and  desired  attributes: 

1.  The  program  must  be  documented  and  tested  In  accordance  with  DoD 
Standard  7935,  dated  15  Feb  1983^,  and  with  the  AFGWC  software  standards'^. 

2.  The  program  must  be  coded  In  ANSI  Standard  X3.9-1978,  as 
Implemented  In  the  UNI  VAC  ASCII  FORTRAN. 
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3.  Any  internal  model  of  Ng{840)  used  should  be  a  continuous  function 
of  latitude,  longitude,  time,  season  (day-of-year) ,  and  sunspot  number  (solar 
epoch) . 

4.  Ideally,  any  internal  Ne(840)  model  should  be  based  on  satellite 
observations.  Provisions  for  updating  this  model  should  be  Included. 

5.  Oata  from  the  SSIE  sensor  should  be  fetched  from  a  user-defined 
section  of  the  Astrogeophysical  Oata  Base  (AGOB). 

6.  Provisions  for  future  incorporation  of  TEC  data,  from  either  the 
AGDB  or  from  a  TEC  data  preprocessor ,  should  be  Included. 

7.  The  analysis  algorithm  should  use  all  data  available  and  be  able  to 
process  data  from  both  northern  and  southern  hemispheres . 

8.  If  no  data  are  available  for  a  particular  program  run,  the  analysis 
should  default  to  the  internal  model. 

9.  The  analysis  should  incorporate  a  "memory"  capability;  i.e.,  the 
analysis  field  from  a  previous  run  could  be  used  as  the  starting  point  for  the 
Current  run. 

10.  The  final  analysis  field  should  not  contain  distinct  features 
(bull's-eyes  or  ridges)  at  points  where  data  are  available. 

11.  The  locations  of  the  grid  points  for  the  output  grids  should  be 
defined  such  that  they  may  be  easily  changed. 

Altnough  this  Is  not  an  officially  approved  statement  of  requirements,  It 
provides  a  starting  point  for  an  investigation  of  alternative  analysis 
methods  for  an  SSIE  Ne(840)  data  preprocessor. 

4.1  Task  4:  Ng(840)  Data  Preprocessor  Improvements 

The  starting  point  for  this  study  was  a  general  review  of  various 
analysis  methods  developed  for  use  In  other  ionospheric  applications,  in 
particular,  those  developed  in  constructing  the  ITS78  model  ionosphere^’^, 
those  used  in  AFGWC  program  UKFIlE  for  updating  the  ITS78  coeff Iclents^.  and 
the  global  analysis  method  used  in  the  4D  model^.  From  this  review,  and  the 
requirements  listed  earlier,  the  following  general  design  was  developed  for 
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an  Improved  SSIELO  program: 

1.  The  core  of  the  program  would  be  a  spectral  model  of  Ne ( 840 )  which 
would  Include  continuous  variability  with  location,  time,  season,  and  solar 

epoch.  The  format  of  this  model  would  be  similar  to  the  ITS78  f  F2  and  M3000 

25  0 

spectral  models  . 

2.  Input  SSIE  Ne( 840)  data  would  be  used  to  modify  the  coefficients  of 
this  model  in  two  stages: 

a.  All  mid-latitude  data  would  be  used  to  adjust  the  model 

through  derivation  of  an  effective  10.7  cm  solar  flux  ( F 10 )  defined  as  the  F 10 

value  which,  when  used  by  the  N^( 840 )  model,  minimizes  the  mean  difference 

between  the  model  and  the  observations.  A  similar  method  for  adjusting  the 

77 

ITS78  f  F2  coefficients  has  been  used  at  AFGWC  for  several  years  . 
o 

b.  All  data  would  be  used  to  modify  the  coefficients  resulting 
from  the  effective  F 10  analysis,  or  from  a  previous  analysis,  using  a  spectral 

OQ 

data  assimilation  method  developed  by  Flattery110  and  used  in  the  global  analy¬ 
sis  section  of  the  4D  model7. 

3.  The  required  UT  grids  would  be  generated  from  either  the  base 
Ne( 840 )  model,  the  results  of  the  effective  F10  nalysis,  the  results  of  the 
global  analysis,  or  from  the  results  of  some  previous  analysis,  at  the  discre¬ 
tion  of  the  user. 

Investigation  of  this  design  would  be  accomplished  in  four  phases: 

1.  Define  the  form  of  the  spectral  Ne(840)  model.  This  would  include 
selection  of  coordinate  system,  functions  to  be  used  in  the  spectral  expan¬ 
sion,  truncation  limits  on  the  spectral  expansion,  and  selection  of  an  initial 
model  of  Ng(840)  on  which  to  base  the  spectral  model.  The  final  step  of  this 
phase  would  be  the  generation  of  coefficients  for  the  spectral  Ng(840)  model. 

2.  Develop  the  effective  F 10  analysis  technique.  The  starting  point 
for  this  phase  would  be  the  AFGWC  UKFIlE  program^. 

3.  Develop  the  global  assimilation  technique.  The  starting  point  for 
this  phase  would  be  the  global  analysis  section  of  the  4D  model7. 

4.  Test  the  full  analysis  program  on  both  model  and  "real-world" 
Inputs. 
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4.1.1  Spectral  Ng(840)  Model 

As  described  in  Scientific  Report  No.  1,  the  functional  form  of  the 
spectral  Ng(840)  model  is 

L  M  N 

fp(\i,t,F,D;840)  =  £  £  ^^(F.D)^  (t)Tm(e  )P^/2(sinA)  (22) 

2.=  1  m=  1  m 

p 

.""(F^D)  =  £  8 lrmn(6 )  Fr  (23) 

<■  r=0  r 

I 

.-;"n(D)  «  £  ‘Hn  Ti(0)  (24) 

i  =0 

where  f  is  the  plasma  frequency  at  840  km  In  MHr;  A  and  <♦,  are  the  F-layer  apex 
latitude  and  longitude;  t  Is  Greenwich  Mean  Time  (UT);  F  is  the  10.7  cm  solar 
radio  flux;  0  is  the  day  of  the  year;  P^2(s1nA)  are  the  associated  Legendre 
functions  of  degree  n  and  order  Sf;  a*™  are  the  model  coefficients;  and  T^(t), 
T(T)(  '),  and  T i ( D )  are  orthonormal  trigonometric  functions  given  by 


T2 j -1  ( x )  =  "=■  cos  Jx’ 

T2j(x)  1  “  Sin  jx. 

The  truncation  limits  for  the  series  In  Equations  22,  23,  and  24  are  L *1 3, 
M=21,  N®20-^,  R=2,  and  1=8.  The  plasma  frequency  from  Equation  22  Is  con¬ 
verted  to  Ne(340)  by 

Ne(840)  =  1.24xl04  f2  el/cm3.  (25) 
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Two  complete  sets  of  model  coefficients,  a  t  have  been  generated  -  one 
using  the  DHR  profile  model  and  one  using  the  Bent  profile  model.  Both  were 
generated  using  the  same  fQF2  and  M3000  values  from  the  ITS78  model  using  the 
procedures  described  in  Section  4.2.1  of  Scientific  Report  No.  1,  with  one 
modification.  It  was  discovered  that  values  of  f  near  the  pole  showed  a  drop 
similar  to  that  found  in  the  SSIELO  model,  even  with  data  values  at  *80 
latitude.  This  was  corrected  by  changing  the  latitudinal  spacing  from  every 
10°  to  a  spacing  defined  by  sin  *  .  =  x^,  i=l,24,  where  the  xi  are  the 
gaussian  quadrature  points  on  the  interval  -1  to  +1.  The  OHR  model  coeffi¬ 
cients  are  used  in  the  analysis  program;  the  Bent  model  coefficients  were 
generated  for  use  in  constructing  data  sets  for  testing  the  analysis  program. 
The  differences  between  these  two  models  is  illustrated  in  Figure  13,  which 
shows  the  northern  geographic  hemisphere  fields  of  Ne(840)  for  a  sunspot 
number  of  120  (F10. 7=164  janskys),  Oay  199  (18  July),  as  generated  from 
Equations  22-24  from  the  DHR  (upper)  and  Bent  (lower)  coefficients.  The  Bent 
coefficient  set  was  generated  when  it  was  discovered  that  test  data  con¬ 
structed  using  the  DHR  coefficients  were  fit  exactly  in  the  effective  F10 
analysis,  leaving  no  residual  error  to  be  fit  in  the  global  assimilation. 


4.1.2.  Effective  F 10  Analysis 

The  objective  of  the  effective  F10  analysis  is  to  find  the  F10  value 
which,  when  used  in  Equation  22,  provides  an  analysis  field  which  has  a  zero 
mean  error  with  respect  to  the  input  data,  i.e.,  Af  =0.  This  value  for  F10  is 
denoted  the  effective  F10,  or  EF10,  for  the  analysis.  The  method  used  is  an 
iterative  one,  the  secant  method,  in  which  subsequent  estimates  of  EF10  are 
calculated  from 


EF 10  =  F.F10 


where  the  superscripts  o  and  -  denote  the  values  from  the  previous  two  itera¬ 
tions  and  Af  is  the  mean  error  between  the  input  data  and  the  field  as 
specified  using  the  EF 10  for  that  iteration.  Once  started,  this  method 
converges  very  rapidly;  however,  it  requires  an  initial  and  one  subsequent 


estimate  for  EF10  to  start. 
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The  Initial  estimate  for  EFIO  is  calculated  by  taking  advantage  of  the 

form  of  the  model  for  f  (840)  given  in  Equation  22.  The  mean  value  for  f  , 

P  P 

averaged  over  all  > ,  \  and  t,  is  given  by 


^P 


(27) 


Assuming  that  the  input  data  are  evenly  distributed  in  t,  ,k,  and  the  mean 
of  the  input  data,  f  ' ,  can  be  used  to  estimate  vj*  by 


Y  =  2*  f  • 

e  P 

From  Equation  23,  yJ*  varies  with  F10  as 

y}1  =  &  +  3.  F10  +  3,  F102 
1  o  1  2 


(28) 


where  the  superscripts  have  been  dropped  on  P.  The  first  estimate  for  EF10  is 
then  defined  as  the  value  of  F10  which  produces  =  ie  using  Equation  28. 
This  leads  to  the  solution 


EF10  =  - 


(29) 


As  the  Input  data  are  not,  as  a  rule,  evenly  spaced,  the  EF10  value  calculated 
from  Equation  29  does  not  necessarily  produce  the  desired  result  of  -fps0. 
Therefore,  further  estimates  are  required. 


The  second  estimate  for  EF10  is  made  using  the  Newton  form  of  Equation 
26.  In  this  method,  the  second  estimate  is  given  by 


in  which  the  derivative  of  Af  with  respect  to  F10  is  required.  Since  the 
Input  fp  values  do  not  vary  with  F10,  an  estimate  for  this  derivative  can  be 
obtained  from  Equations  27  and  28  yielding 


-  -if 

?FT0  “  0 


F10 


/2 


:1  *  2F1°  *2) 
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so  Equation  30  becomes 


EF10+  =  EF10° 


l-Jl  Af° 

_ He _ 

+2  EF 1 0°  S2)  * 


(31) 


As  with  Equation  29,  this  will  provide  only  a  next-best  estimate  for  EF10,  as 
it  also  uses  the  assumption  that  the  data  are  uniformly  distributed  in  space 
and  time,  and  further  iteration  using  Equation  26  may  be  required.  In  actual 
testing  (see  Section  4.1.5),  the  second  estimate  was  within  the  stop  criterion 
of  _'.f  -  0.005MHz  in  9  of  12  cases,  and  only  one  further  iteration  was  required 
in  the  remaining  3  cases. 

The  effective  F10  analysis  may  also  be  run  using  a  set  of  coeffi¬ 
cients  produced  by  a  previous  analysis.  Equation  24  can  be  rewritten 


i‘mn  =  mn  _  |,imn  pl0  +  g|nin  F1Q2j  _ 

By  using  the  and  EF10  values  from  a  previous  analysis  in  this  equation, 
the  ^omn  coefficients  are  reset,  and  the  effective  F10  analysis  can  then  be 
done  using  the  updated  set  of  S*"1"  coefficients. 


4.1.3.  Global  Analys i s 


The  global  analysis  technique  is  identical  to  that  employed  by  the  AWS  4D 
ionospheric  analysis  model**.  This  technique,  developed  initially  for  mete- 
orological  applications  ,  modifies  the  coefficients  of  a  spectral  repre¬ 
sentation  of  a  field  such  that  the  new  coefficients  produce  a  field  that  fits 
input  data  where  data  are  available  and  merges  smoothly  into  the  original 
field  away  from  the  data  points.  In  this  analysis,  the  >™n  coefficients  of 
Equation  22  are  iteratively  updated  using 


mn ' 


mn 


w£V 

j  =  l 


mn 


<y 


V 


tVC0S,j 


(32) 


A*n( 


m/2, 


rtV  *  TTtj>V:J)pn 
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where  J  Is  the  number  of  data  points,  A"f p(  A  '"j * j ^  *s  the  di^^erence  between 
the  Input  f  and  f  calculated  using  the  >  ^  coefficients  in  Equation  22  at 
the  jth  data  point,  and  W  is  a  weight  factor  given  by 

W  =  w^  A  A AC A  t 

where  wc  is  a  convergence  control  parameter  and  A'-,  A:,  :.t  define  the  lati¬ 
tude,  longitude,  and  time  spacing  of  the  fictitious  grid  over  which  the 
analysis  is  defined®.  The  Iteration  continues,  until  a  desired  RMS  difference 
between  the  data  and  the  analysis  RMS  difference  and  begins  to  diverge. 

In  order  to  control  the  introduction  of  spurious  small-scale  features 
into  the  analysis,  the  analysis  is  designed  to  begin  the  iterative  update  on 
the  lower-order  coefficients  only.  This  is  done  by  using  values  for  L  and  M 
in  Equation  22  lower  than  the  maximum  used  in  the  model  (currently  13  and  21), 
and  increasing  these  values  with  each  Iteration  until  the  maximum  truncation 
limits  are  reached.  By  starting  the  analysis  without  the  higher-order  terms 
in  me  expansion,  the  difference  between  the  data  and  the  initial  field  will 
be  reduced  first  through  modification  of  the  larger  scale  features  of  the 
field,  and  only  the  residual  differences  will  go  into  smal ler-scale  features. 
Initially,  it  was  planned  to  reduce  only  the  longitude  (M)  and  t;me  (L) 
truncation  limits,  and  this  is  what  was  used  in  the  developmental  testing 
described  in  Section  4.1.4.  However,  during  the  study  of  using  No(840)  values 
derived  from  TEC  dati,  described  in  Section  4.2,  it  was  decided  to  allow  the 
latitude  (N)  truncation  limit  such  that  the  value  of  N  at  m=l  was  N=M-1,  thus 
preserving  the  modified  triangular  truncation  scheme. 

4.1.4  Preprocessor  Program  (GRIONE)  Description 

The  analysis  procedures  described  in  Sections  4.1.2  and  4.1.3  were 
Implemented  in  a  test  version  of  an  Ng( 840)  data  preprocessor  (program  GRIDNE) 
designed  and  coded  to  meet  the  requirements  worked  out  with  AFGWC/SDDE  and 
AFGl/PHG.  Figure  14  is  an  organization  chart  for  program  GRIDNE  which 
illustrates  the  Internal  linkages  and  relationships  between  the  various 
routines  which  make  up  the  program.  A  short  description  of  the  function  of 
each  of  the  routines  listed  in  Figure  14  is  presented  in  Table  4,  and  a 
complete  listing  of  all  non-AFGWC  routines  is  included  in  Appendix  C. 


Figure  14.  Program  GRIDNE  Structure  Chart, 


Table  4.  Description  of  program  GRIDNE  routines. 


ADJF10  -  Calculate  the  next  estimate  of  the  effective  F10  using  the  secant 
iteration  method  (Implements  Equation  26). 

‘AGFECH  -  Fetch  SSIE  data  records  from  the  AFGWC  Astrogeophysical  Data  B$se 
(AGDB) . 

*AGDATE  -  Determine  the  active  dates  and  julian  hours  on  the  AGDB. 

AGDCHK  -  Determine  if  an  input  date  and  time  is  active  on  the  AGDB. 

ALFMTT  -  Calculate  the  associated  Legendre  functions  for  a  modified  triangle 
truncation  scheme. 

ANALNE  -  Perform  the  global  analysis  (implements  Equation  32). 

APxTAB  -  Calculate  the  F-layer  apex  coordinates  of  an  input  geographic  lati¬ 
tude  and  longitude. 


Table  4.  (Continued) 


CALCNE  -  Calculate  a  value  for  f  (840)  and  N  (840)  for  a  given  time  and 
location  (implements  Equation  22).  e 

C0EF10  -  Input  or  output  a  set  of  y™"  coefficients  for  Equation  22. 

OTGCHK  -  Error  check  a  user- input  date  and  time. 

GETCO  -  Input  the  coefficients  for  Equation  24  and  sum  over  the  date 

( implements  equation  24). 

GETNE  -  Fetch  SSIE  data  records  from  the  AGDB,  erro  check  the  records,  and 
store  them  on  a  temporary  file. 

♦GRIDIO  -  Write  Ng(840)  grid  sets  to  file  FOURDGRIDS. 

GRIDNE  -  Main  driver  routine. 

IN  I T  -  Initialize  all  variables  requiring  initialization  and  load  the  array 
which  defines  the  FOURDGRIDS  output  grid. 

INPUT  -  Interact  with  the  user  for  all  required  user  inputs,  fetch  SSIE  data 
from  the  AGDB  and  prepare  the  data  for  analysis  if  necessary,  and 
fetch  the  necessary  coefficient  sets. 

INTERP  -  Interpolate  in  the  F-1ayer  apex  coordinate  look-up  table. 

OFFSET  -  Calculate  the  offset  dipole  coordinates  of  a  geographic  location. 

0UT6RD  -  Build  and  output  Ng(840)  grids  to  file  FOURDGRIDS. 

SETFIO  -  Perform  the  effective  F10  analysis  (implements  Equation  29). 

SUMF10  -  Sum  the  2^  coefficients  over  F 10  (implements  Equation  23). 

TFUNC  -  Calculate  the  orthononnal  trignometric  series. 

THPIN  -  Input  SSIE  Ne(840)  records  from  a  temporary  file. 

TMPOUT  -  Output  SSIE  N  (840)  records  to  a  temporary  file. 


YESNO 


Query  the  user  for  a  yes/no  answer. 


♦AFGWC  Routines 


I 


In  tests  at  AFGWC  on  a  SPERRY  11/82  computer,  the  program  size  was 
determined  to  be  roughly  55K  words,  and  the  CPU  time  required  for  an  analysis 
run  can  be  estimated  by  the  algorithms 

CPUj  =  3.8  np  ♦  12.7  ng  -  0.5  seconds 

for  processing  data  from  one  satellite  for  a  24-hour  period  (1351  data 
points)  and 

CPU^  =  7.1  np  ♦  23.8  n^  -  8.0  seconds 

for  data  from  two  satellites,  where  np  Is  the  number  of  Iterations  required  to 
determine  the  effective  FI 0  value  and  n.  is  the  number  of  Iteration.,  In  the 
global  analysis.  From  the  tests  conducted,  it  appears  that  nominal  values  for 
np  and  nQ  are  2  and  6,  respectively,  for  a  two-satellite  CPU  time  of  roughly 
150  seconds.  Worst-case  values  were  np  =  3  and  n^  *  10  for  a  two-satellalte 
CPU  time  of  roughly  250  seconds. 

Inputs  from  the  user  are  as  follows: 

a.  Oate  ( YYMMDD)  and  sunspot  number  for  the  analysis. 

1  b.  The  maximum  number  of  Iterations  for  the  effective  F10  analysis, 

and  the  value  of  the  average  Afp  at  which  to  stop  the  analysis. 

c.  The  maximum  number  of  iterations  for  the  global  analysis;  the  value 

'  of  the  RMS  if  at  which  to  stop  the  analysis;  the  initial  truncation  limits  to 

be  used  for  L  and  M  in  Equations  22  and  32;  the  spacing  of  the  latitude, 
longitude,  time  grid  over  which  Equation  32  Is  defined;  and  the  value  for  the 
j  convergence  control  parameter,  Wc, 

d.  Various  program  control  parameters  specifying  the  source  of  the 
data  to  be  analyzed,  the  coefficients  to  use  In  defining  the  Initial  analysis 
field,  which  analyses  to  run,  and  what  items  to  output. 

!  The  program  uses  two  disk  storage  files,  NE840FILE  and  NETEMPFILE.  File 

j  NE840FILE,  roughly  103K-words  in  size,  contains  the  F-layer  apex  coordinate 

|  look-up  table  (8190  words),  the  coefficients  which  define  the  Ne(840) 

!  model  (77,571  words),  five  sets  of  y?n  coefficients  from  previous  analyses 

I  (2873  words/set,  14,365  words  total),  and  six  28-word  file  control  sectors, 

j  Subroutine  APXTA8  reads  in  the  apex  coordinate  table,  subroutine  GETCO  reads 

J 

i 

t 

i 

i 


i'i+f v 


f 
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In  the  coefficients,  and  subroutine  C0EFI0  reads  and  writes  the  '  T n 

coefficients  sets.  File  NETEMPFILE  Is  a  temporary  file  used  for  external 
storage  of  the  Ng(840)  data  to  be  analyzed.  Subroutine  GETNE  strips  the 
Ne(840)  data  from  SSIE  AGOB  records  and  builds  this  file  via  subroutine 
TMPOUT.  Subroutines  SETFljJ  and  ANALNE  read  the  Ne(840)  data  from  the  file  via 
function  TMPIN. 

4.1.5  Preprocessor  Program  Tests 

Tests  designed  to  determine  the  capabilities  and  limitations  of  the 
GRIDNE  program  were  run  at  AFGWC  during  the  period  10-14  September  1984.  A 
total  of  three  tests  were  conducted: 

a.  Test  1  -  Test  the  overall  capabilities  and  limitations  of  the 
analysis  methods  for  a  range  of  sunspot  numbers. 

b.  Test  2  -  Test  the  effects  of  varying  the  convergence  control  para¬ 
meter,  Wc . 

c.  Test  3  -  Test  the  effects  of  using  various  initial  truncation  limits 
for  M  and  L  in  Equations  22  and  32. 

The  input  parameters  used  in  the  tests  are  listed  in  Table  5. 

4. 1.5.1  Description  of  Test  Data  and  Procedures 

All  tests  were  run  using  N^( 840)  data  generated  by  using  the  Bent 
coefficients  in  Equations  22-25.  Two  sets  of  data  were  generated  for  each 
test,  one  to  be  used  as  Input  to  program  GRIDNE  and  one  to  be  used  to 
determine  how  well  the  GRIDNE  analyses  performed.  The  first  set,  denoted  the 
satellite  data  set,  consists  of  simulated  SSIE  Ng(840)  data  for  a  24-hour 
period  (1351  data  points)  for  one  or  two  DMSP  satellites  with  ascending  modes 
at  0600LT  and  10001T.  The  second  set,  denoted  the  Bent  grid  set,  consists  of 
24  Ne(840)  grids,  one  grid  per  UT  hour,  covering  the  northern  geographic 
hemisphere  from  the  equator  to  80°N  latitude  in  10°  increments  and  03  to  345°E 
longitude  In  15°  increments.  A  test-data  generating  program,  ORBTNE,  was 
developed  which  can  generate  either  data  set  or  both  for  a  given  date  and 
sunspot  number  (SSN). 

The  test  procedures  were  as  follows: 

a.  Run  test  program  0R3TNE  to  generate  the  satellite  and  grid  data  sets 
for  the  data  and  SSN  specified  in  Table  5. 
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Table  5.  Test  Parameters 


I .  Test  1  -  Analysis  Capabilities 
Oate:  840911 

Initial  truncation  limits:  M»13,  1*7 

Grid  spacing:  10°  latitude,  15°  longitude,  l*1  time 

Convergence  control :  *c*1.0 

Case  la  -  SSN-10,  0600LT  satellite 

Case  lb  -  SSN-10,  0600LT  and  1000LT  satellites 

Case  lc  -  SSN»60,  0600LT  satellite 

Case  Id  -  SSN*60,  0600LT  and  10001T  satellites 

Case  le  -  SSN*120,  0600LT  satellite 

Case  If  -  SSN-120,  0600LT  and  1000LT  satellites 

I I .  Test  2  -  Converqei.ee  Control  Parameter  Effects 
Oate:  840911 

Initial  truncation  limits:  M*13,  1*7 

Grid  spacing:  10°  latitude,  15°  longitude,  lh  time 

Data  set:  SSN-60,  0600LT  and  10001T  satellites 

Case  2a  •  wc*1.50 

Case  2b  -  wc*1.25 

Case  2c  -  *c*1.10 

Case  2d  -  *c*0.90 

Case  2e  -  *c*0.80 

Case  2f  -  w  =0.50 
c 

III.  Test  3  -  Truncation  Limit  Effects 
Date:  840911 

Grid  spacing:  10’  latitude,  15’  longitude,  time 
Convergence  control:  1.0 

Data  set:  SSN=60,  0600LT  and  1000LT  satellites 

Case  3a  -  M-13,  1=5 

Case  3b  -  M=15,  1*7 

Case  3c  -  M=17,  1=9 

Case  3d  -  M*19,  1=11 

Case  3e  -  M°21,  L*13 
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b.  Run  program  6RIDNE  to  process  the  satellite  data  set.  Outputs 
Include  goodness-of-f It  statistics  from  both  the  effective  F10  and  global 
analyses,  and  the  v1™  coefficients  from  the  Initial  analysis  field  (calcu¬ 
lated  from  Equations  23  and  24  using  the  PHR  coefficients  and  the  same 
date  and  SSN  used  by  program  ORBTNE),  the  effective  F 10  field  (calculated  from 
Equation  23  using  the  effective  F10  value),  and  the  global  analysis  field 
(calculated  Iteratively  from  Equation  32). 

c.  Generate  a  grid  set  from  each  of  the  three  >  ;mn  sets  output  from 
program  GRIONE.  Compare  each  of  these  three  to  the  Bent  grid  set  and  generate 
the  following  comparison  statistics: 

(1)  A  grid  set  which  contains  the  difference  between  the  analysis 
and  Bent  grid  sets,  denoted  the  ANg(840)  grid  set. 

(2)  The  RMS  ANg(840)  and  RMS  percent  ANg(840)  for  all  24  grids 
(5184  data  points),  and  as  functions  of  latitude,  longitude,  UT,  and  local 
time. 

Figures  15  and  16  are  samples  of  the  analysis  results  which  will  be 
presented  for  the  various  tests.  This  particular  test  was  for  18  July  with  a 
sunspot  number  of  120,  processing  data  from  two  DMSP  satellites.  In  Figure  15, 
the  upper  plot  shows  the  variation  of  RMS  -Ng(840)  with  latitude  for  the 
Initial  analysis  field  (dotted  line),  the  effective  F10  (EF10)  field  (dashed 
line),  and  the  global  analysis  field  (solid  line).  The  lower  plot  In  this 
figure  shows  the  variation  of  ANg(840)  with  local  time  for  all  three  fields. 
The  vertical  dotted  lines  In  this  figure  Indicate  the  local  times  of  the 
ascending  (0600  and  1000)  and  descending  (1800  and  2200)  nodes  of  the  DMSP 
orbits  used  In  the  analysis. 

Figure  16  shows  contour  plots  of  the  0000UT  analysts,  Ne(840),  and  Bent 

grids.  The  three  pairs  of  plots  are,  from  the  top  down,  tne  Initial  analysis 

grid  and  Its  ^(840)  grid,  the  EF10  analysis  grid  and  its  ANe(840)  grid,  and 

the  global  analysis  grid  and  its  ANg(840)  grid.  The  single  plot  at  the  bottom 

Is  the  Bent  grid.  Ideally,  If  the  analysis  program  worked  perfectly,  the 

global  analysis  grid  and  the  Bent  grid  would  be  identical.  (Note:  The 

decision  to  display  the  0000UT  grids  as  representative  grids  for  the  entire 

24-hour  grid  set  was  Initially  arbitrary,  but  In  reviewing  the  analysis 

results  It  turns  out  that  the  total  RMS  ‘N  (840)  for  the  0000UT  grid  was  close 

to  the  RMS  AN  (840)  for  the  entire  set  and  never  represented  either  the  best- 
e 

or  worst-case  geld. ) 
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4. 1.5.2  Test  Results 

The  analysis  goodness-of-f it  statistics  from  all  Test  l  test  cases  are 
suirmarized  in  Table  6  (EF10  analysis)  and  Figure  17  (global  analysis).  Table 
6  lists  the  initial  EF10  estimate,  calculated  from  Equation  29;  the  second 
EF10  estimate,  calculated  from  Equation  31;  and  the  EF10  estimate,  calculated 
after  one  iteration  of  Equation  26.  Included  in  this  table  Is  the  average 
difference  between  the  analysis  field  and  the  input  data,  £fp,  for  each  EF10. 
In  all  cases,  the  EF10  estimated  after  one  iteration  of  Equation  26  produced  a 
zero  mean  difference,  as  desired.  As  Figure  17  shows,  the  global  analyses  all 
converged  uniformly,  reaching  a  point  of  diminishing  returns  at  roughly  itera¬ 
tion  number  6  or  7. 


Table  6.  Results  of  the  effective  F10  analyses  in  Test  1. 


Equation  29  Equation  31  First  Iteration 


Test 

Case 

EF10 

EF10 

EF10 

la 

69.78 

0.0471 

66.87 

0.0043 

66.58 

0 

lb 

74.26 

0.0501 

71.16 

-0.0010 

71.22 

0 

lc 

90.46 

0.0160 

89.45 

0.0013 

89.36 

0 

Id 

97.46 

0.0568 

93.87 

-0.0013 

93.95 

0 

le 

119.31 

-0.0205 

120.63 

0.0014 

120.73 

0 

If 

130.09 

0.0664 

125.77 

0.0014 

125.86 

0 

Figure  18  illustrates  how  well  the  input  data  are  fit.  The  two  plots  In 
this  figure  show  the  f p ( 840)  data  input  for  the  first  full  orbit  from  both 
satellites  (solid  lines)  from  test  case  If  (SSN=120,  two  satellites)  and 
fp(840)  calculated  for  the  same  orbit  locations  from  the  initial,  EF10,  and 
global  analysis  fields  (dotted  lines).  The  final  RMS  error  for  the  entire 
analysis  was  0.060  MHz  in  plasma  frequency  and  3.8x10^  el/cm^  In  electron 
density.  This  level  of  accuracy  is  reflected  in  the  small  deviations  of  the 
global  analysis  curves  from  the  input  curves  in  Figure  18.  The  only  major 
deviations,  occurring  near  orbit  data  point  numbers  28  and  75,  are  near  the 
F-layer  apex  poles.  This  is  a  difficult  region  to  make  changes  to,  as  all 
data  are  weighted  by  the  cosine  of  the  F-layer  apex  latitude  in  the  global 
analysis  (see  Equation  31),  causing  data  near  the  poles  to  be  almost  com¬ 
pletely  ignored. 


54 
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0600LT  SATELLITE  -  ORBIT  1 


ORBIT  DATA  POINT  NUMBER 


1000LT  SATELLITE  -  ORBIT  1 


ORBIT  DATA  POINT  NUMBER 


Figure  18.  Analysis  results  for  the  first  full  orbit  of  both  DMSP 
satellites  for  test  case  If. 


Table  7  provides  a  summary  of  the  RMS  ANe(840)  and  RMS  percent  ANg(840) 


Table  7.  Sumnary  of  the  results  of  the  Ng(840)  grid  set  comparisons. 

RMS  Ne(840)  RMS  X  Ne(840) 

Test  ~  ’  '  .  *  ~  ~  ” 


Case 

Initial 

EF10 

Global 

Initial 

EF10 

Global 

la 

0.99 

1.18 

0.89 

58.0 

61.1 

53.2 

lb 

0.99 

1.29 

0.77 

58.0 

60.8 

41.2 

lc 

2.82 

2.08 

1.76 

68.3 

59.1 

51.7 

Id 

2.82 

2.14 

1.40 

68.3 

58.6 

39.2 

le 

5.77 

3.63 

3.30 

78.2 

55.4 

49.9 

If 

5.77 

3.59 

2.58 

78.2 

54.7 

39.3 

xlO4  el/cm^ 

Percent 

statistics  for 

the  N  (840) 

grid  set  comparisons  for 

all  six 

test  cases. 

In 

each 

set  of  three  numbers. 

the  first 

is  from  comparing  the 

initial  field 

to 

the 

Bent  field, 

,  the  second 

from  the  1 

EF10  field  and 

the  Ben 

t  field,  and 

the 

third  from  the 

global  analysis  field 

to  the  Bent 

field. 

In  general. 

the 

following  observations  can  be  made  from  these  comparisons: 

a.  Not  surprisingly,  the  2-satellite  cases  (lb.  Id,  If)  provided 
better  overall  fits  to  the  Bent  grid  sets  than  the  corresponding  1-satellite 
cases  (la,  lc,  le).  The  1-satellite  cases  tend  toward  a  final  RMS  percent 
ANe(840)  of  around  50X,  while  the  2-satellite  cases  tend  toward  40X. 

b.  In  all  cases,  the  final  RMS  ANg(840)  and  RMS  percent  ANg(840)  are 
lower  than  that  obtained  with  the  Initial  analysis  field,  and  the  total 
processing  gain  as  reflected  in  both  figures  increases  with  sunspot  number. 

c.  The  effectiveness  of  the  EF10  analysis  increases  dramatically  with 
sunspot  number.  There  Is  actually  a  processing  loss  for  low  sunspot  numbers, 
at  least  In  terms  of  RMS  differences,  in  the  EF10  analysis.  This  is  largely 
due  to  the  fact  that  the  EF10  analysis  tries  to  reduce  the  average  difference 
to  zero,  which  can.  In  certain  circumstances.  Increase  the  RMS  difference.  The 
fact  that  the  analysis  is  done  in  plasma  frequency  wnile  the  statistics  in 
Table  7  are  in  electron  density  can  also  contribute  to  this  effect. 
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The  observations  are,  for  the  most  part,  supported  by  the  more  detailed 
analysis  results  shown  in  Figures  19-30.  Two  figures  are  presented  for  each 
of  the  six  test  cases  -  the  plots  of  RMS  ANe(840)  as  a  function  of  latitude  and 
local  time  (as  in  Figure  15)  first,  followed  by  the  OOOOUT  grid  contour  plots. 
The  main  observation  that  can  be  drawn  from  these  figures  is  that  the  global 
analysis  field  reverts  back  to  the  EF10  analysis  field  within  roughly  30°  of 
longitude,  or  2h  of  local  time,  of  an  orbit  track.  In  other  words,  It  Is  the 
relative  accuracy  of  the  shape  of  the  Ng(840)  model  used  initially  to  generate 

the  coefficients  used  In  Equation  24  (the  DHR  model  in  the  current  tests) 

ri 

and  the  absolute  accuracy  provided  by  the  EF10  analysis  which  determines  the 
accuracy  of  the  field  more  than  30°  or  2h  away  from  data  points.  Conversely, 
where  there  is  an  adequate  collection  of  data  the  global  analysis  can  repro¬ 
duce  the  actual  field  quite  accurately,  almost  in  spite  of  the  field  that  the 
global  analysis  starts  with.  This  can  be  seen  In  the  dramatic  Improvement  In 
the  RMS  lNe(840)  versus  local  time  plots  at  around  16-20LT. 

The  objective  of  Test  2  was  to  study  the  effects  of  varying  the  conver¬ 
gence  control  parameter,  Wc,  on  the  global  analysis  algorithm  (Equation  32). 
The  $$N=60,  2-satellite  data  base  from  test  case  Id  was  used  for  this  test. 
The  results  are  summarized  in  Table  8,  which  lists  the  RMS  Afp(840)  values  for 
each  iteration  step  and  the  RMS  ^Ne(840)  and  RMS  percent  ANg(840)  values  from 
the  Ng(840)  grid  set  comparisons.  These  results  are  somewhat  surprising  in 
that  the  relatively  arbitrary  choice  of  A*  =  10°  ,  Ao=  15°,  ~>t  a  l^1  and  Wc=1.0; 
values  which  were  carry-overs  from  previous  work  on  the  40  analysis  model; 
provide  the  best  fit  in  the  Ng(840)  grid  comparisons.  A  limited  check  of 
these  results  using  the  SSN-60,  1-satellite  data  base  confirmed  the  finding 
that  this  choice  of  values  for  the  parameters  which  go  into  calculating  the 
value  for  W  in  Equation  32  appears  to  be  close  to  optimal  for  processing 
orbital  data.  (In  running  the  tests  described  in  section  4.2,  In  which 
roughly  5000  data  points  were  processed,  Wc  had  to  be  reduced  to  roughly  0.6 
to  obtain  convergence.) 
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RMS  ERROR  NEC840?  RMS  ERROR  NE<840> 


X104  RMS  ANALYSIS  ERROR  SUMMARY 


LATITUDE 

X104  RMS  ANALYSIS  ERROR  SUMMARY 


LOCAL  TIME 


Figure  19.  Summary  of  the  Ne(840)  grid  set  coinoari sons  for  test  case 
la  ( SSN  =  10,  1  satellite). 
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RMS  ERROR  NEC840:>  RMS  ERROR  NEC840) 


LATITUDE 


LOCAL  TIME 

gure  21.  Su  unary  of  the  Ne(840)  grid  set  comparisons  for  test  case 
lb  ( SSN  =  10,  2  satellites). 


RMS  ERROR  NEC840)  RMS  ERROR  NEC840) 


LATITUDE 


LOCAL  TIME 


Figure  23.  Summary  of  the  Ne(840)  grid  set  comparisons  for  test  case 
lc  ( SSN=60,  1  satellite). 
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Figure  24.  Analysis  results  for  the  0000UT 
grid  for  test  case  1c  (SSN=60, 

1  satellite). 
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Table  8.  Test  2  results  summary. 


I 

T 

E 

R 

A 

r 

i 

o 

N 


2a 

2b 

2c 

Id 

2d 

2e 

2f 

'»se 

1.50 

1.25 

1.10 

1.00 

0.90 

0.80 

0.50 

r 

1  0.577 

0.577 

0.577 

0.577 

0.577 

0.577 

0.577 

2  0.352 

C.311 

0.332 

0.339 

0.350 

0.366 

0.431 

3  0.412 

0.259 

0.219 

0.214 

0.224 

0.243 

0.330 

4  -0- 

0.269 

0.156 

0.127 

0.128 

0.144 

0.235 

5 

-D- 

0.139 

0.087 

0.079 

0.089 

0.161 

6 

0.137 

0.069 

0.058 

0.066 

0.120 

MHz 

7 

0.140 

0.058 

0.046 

0.052 

0.095 

8 

-0- 

0.049 

0.038 

0.043 

0.079 

9 

0.043 

0.032 

0.037 

0.067 

10 

0.039 

0.028 

0.032 

0.058 

RMS  - 

1.40 

1.42 

1.43 

1.51 

X 13* el /Cm3 

RMSX  - 

- 

- 

39.2 

39.8 

40.1 

41.2 

t 

The  objective  of  Test  3  was  to  study  the  effects  of  varying  the  initial 
truncation  limits  on  the  longitude  (M)  and  time  (L)  expansions  in  Equations  22 
and  32.  As  in  Test  2,  the  SSN=60,  2-satellite  data  base  from  test  case  Id  was 
used  for  this  test.  The  results  of  this  test  are  summarized  in  Table  9,  which 
lists  the  RMS  — f p ( 840 )  values  for  each  iteration  step  and  the  RMS  .'.Ne(340i  and 
RMS  percent  ^g(840)  values  from  the  Ng(840)  grid  set  comparisons.  The  gain 
achieved  by  starting  the  analysis  with  only  the  lower  wave  numbers  is  very 
obvious  in  these  results.  Although  the  cases  using  higher  truncation  limits 
initially  do  better  at  fitting  the  input  data,  all  cases  are  roughly  equal  by 
iteration  step  7,  and  the  lower  truncation  limit  cases  produce  better  fits  at 
iteration  step  10.  Even  more  significant  is  the  gain  evident  in  the  fits  to 
the  overall  fi*ld  reflected  in  the  RMS  Ne(840)  and  RMS  percent  N e 1 G - ? ' 
parameters,  indicating  that  the  information  in  the  input  data  is  being  'nor? 
effectively  spread  away  from  the  data  point  locations  by  using  lowe’*  imf;’ 
truncation  values. 
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Table  9.  Test  3  results  summary. 


3a 

3b 

3c 

3d 

3e 

Case 

1 

i 

2 

T 

13/5 

0.573 

15/7 

0.577 

17/9 

0.580 

19/11 

0.579 

21/13 

0.578 

M/L 

0.359 

0.350 

0.323 

0.259 

0.224 

i 

3 

r 

0.245 

0.223 

0.170 

0.136 

0.136 

4 

r, 

0.163 

0.127 

0.099 

0.093 

0.096 

K 

5 

A 

0.099 

0.079 

0.072 

0.071 

0.073 

MHz 

M 

6 

T  7 
'  3 

0.064 

0.058 

0.055 

0.058 

0.058 

0.049 

0.046 

0.045 

0.050 

0.049 

0.039 

0.038 

0.038 

0.046 

0.04? 

u 

9 

N 

10 

0.033 

0.032 

0.033 

0.046 

0.038 

0.028 

0.028 

0.029 

0.046 

0.036 

RMS 

1.35 

1.42 

1.53 

1.61 

1.77 

xl04el/cm^ 

RMS* 

38.5 

39.9 

41.4 

44.1 

47.3 

% 

An  additional  gain  is  in  processing  time.  Since  there  are  fewer  terms  to 
s um  over  in  both  Equation  22  and  32,  the  lower  the  initial  truncation  limits 
the  faster  the  global  analysis  will  run.  The  CPU  times  for  the  two  extreme 
cases  shown  in  Table  9  were  203  seconds  for  case  3a  and  225  seconds  for  case 
3e,  a  decrease  in  resource  usage  of  10*  in  order  to  achieve  a  better  analysis. 

4.1.6  Limitations 

The  major  weakness  of  the  preprocessor,  as  implemented  for  testing,  is  in 
the  spectral  Ne(840)  model  which  serves  as  the  first-guess  field.  Both  models 
used,  the  AWS  DHR  model  and  the  Bent  model,  were  designed  primarily  for  the 
calculation  of  integrated  electron  content  values,  even  though  the  Bent  model 
was  initially  derived  from  topside  lonosonde  data.  The  effects  of  this  are 
most  apparent  in  the  fields  generated  from  the  DHR  model  coefficients  which 
often  include  odd  features,  such  as  the  two  features  at  105°E  and  240°E  in  the 
upper  panel  of  Figure  13,  induced  by  the  local-time  arying  topside  scale 
height  factors  used  by  the  Dr'X  model.  These  factors  also  produce  the  Intense 
second  diurnal  maximum  at  around  17001T  in  the  equatorial  region.  On 


the  other  hand,  the  Ne(840)  values  calculated  front  the  Bent  coefficients 
appear  to  be  too  low  for  a  given  F10  value.  This  is  less  of  a  problem,  given 
the  results  of  the  effective  F10  tests,  but  it  is  not  clear  that  the  relative 
field  shape  (in  latitude  and  longitude)  provided  by  the  Bent  model  is  any 
better  or  worse  than  that  provided  by  the  DHR  model.  Additionally,  neither 
model  does  well  at  high  latitudes  or  during  geomagnetic  storms. 

A  related  limitation,  also  due  to  limitations  placed  on  the  imple¬ 
mentation  of  the  algorithms  rather  than  the  algorithms  themselves,  is  the 
maximum  truncation  limits  used  in  Eguations  22-24  and  32,  which  limit  the 
inherent  resolution  of  the  analysis.  These  were  selected  by  a  trade-eff 
between  desired  resolution  (in  latitude,  longitude,  and  time)  and  program 
size  and  run-time,  and  should  be  redefined  based  on  actual  user  requirements 
and  resources  in  any  operational  implementation  of  this  preprocessor. 

4.2  Task  5:  Combined  TEC  and  Ne(340)  Preprocessor 

Tests  were  made  using  the  GRIDNE  program  described  in  Section  4.1  to 
investigate  the  effects  of  incorporating  Ng(840)  data  inferred  from  "EC 
observations  into  analyses  of  SSIE  Ne(840)  data.  Three  tests  were  rUn  usinq 
data  generated  using  the  Bent  coefficients  with  Equations  22-24,  and  one  test 
was  run  using  TEC  data  from  the  AFGWC  TEC  data  preprocessor  and  SSIE  data  from 
the  AGOB. 

In  the  first  three  tests,  it  is  assumed  that  a  processor  is  available 
which  will  produce  estimates  of  Ne(840)  from  TEC,  so  that  only  Ng(840)  data  is 
provided  to  the  program.  These  three  tests  use  the  SSN=60,  2-satellite  data 
base  from  the  testing  described  in  Section  4. 1.5. 2,  and  add  three  different 
TEC-inferred  ( 840 )  data  sets: 

1.  Data  from  the  12  AWS  polarimeter  locations  (stations'  listed  in 
Table  10  for  each  of  24  hours  (288  additional  data  points). 

2.  Data  on  an  F-layer  apex  qnd  (Grid  1)  from  2(f  N  to  6C  N  latitude  ;n 
10°  increments,  0°  to  345°E  longitude  in  153  increments,  and  OO'JT  to  ’o 
one-hour  increments  ^2880  additional  data  points). 


73 


Table  10.  TEC  station  penetration  point  locations. 


Station 

Latitude 

Longitude 

Athens 

34°N 

18°E 

Goose  Bay 

49°N 

294°E 

Osan 

35°N 

128°E 

Patrick  AFB 

27°N 

2  78°E 

Ramey 

17°N 

291°E 

Sagamore  Hill 

39°N 

288°E 

Boulder 

37°N 

254°E 

La  Posta 

244°E 

Taiwan 

23°N 

122°E 

Palehua 

20°N 

203°E 

Shemya 

473N 

164°E 

Anchorage 

54  °N 

214°E 

Table  11.  TEC  tests  results  summary. 


EF10 

fp(840) 

RMS 

RMSX 

No  TEC 

93.95 

0.039 

1.40 

39.2 

Stations 

96.56 

0.038 

1.36 

37.3 

Grid  1 

105.93 

0.032 

1.43 

33.6 

Grid  2 

99.92 

0.030 

0.81 

23.3 

Stations  (NT) 

96.56 

0.036 

1.31 

36.2 

Grid  3  (NT) 

101.32 

0.034 

0.78 

20.8 
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3.  Data  on  an  F-1ayer  apex  grid  (Grid  2)  from  0°  to  70~N  latitude  in 
10°  increments  with  the  same  longitude  and  time  coverage  as  in  (2)  (4608  addi¬ 
tional  data  points) . 

The  last  two  data  sets  were  assumed  to  be  from  some  TEC  data  preprocessor . 

The  results  of  these  three  tests  are  summarized  in  Table  11,  and  the 

OOOOirr  Ne(840)  and  *Ne(840)  grids  for  all  three  are  shown  in  Figure  31  along 

with  the  Bent  Ng(840)  grid.  As  Table  3  shows,  the  quality  of  the  analysis 

Improves  with  increasing  data  input,  with  the  exception  of  the  Grid  1  case 

where  the  RMS  -Ne(84C)  increased  although  the  RMS  percent  J.Ne(840)  decreased. 

This  is  due  to  the  higher  EF10  value  calculated  for  this  case  (105.93),  as  can 

be  seen  in  Figure  31.  Although  the  mio  latitudes  are  now  very  close  to  the 

Bent  field  (the  results  of  the  global  analysis),  the  low  latitudes  are  much 

higher  than  in  the  other  two  cases,  a  direct  result  of  the  higher  EF 10  value. 

Since  the  Ne(840)  values  are  higher  at  low  latitudes  than  mid  latitudes,  this 

caused  the  RMS  .tNe(340)  value  to  increase  even  though  the  overall  fit,  as 

reflected  in  the  RMS  percent  N  (840)  value,  was  better.  In  the  Grid  2  case, 

0 

where  almost  an  entire  hemispheric  grid  was  processed,  the  global  analysis  has 
modified  the  field  completely  so  that  it  is  almost  identical  to  the  Bent  grid. 

As  stated  in  Section  4.1.3,  the  variable  truncation  limit  scheme  was 
applied  only  to  the  longitude  and  time  dimensions.  The  latitude  dimension  was 
not  Included,  as  the  orbit-track  data  set  provided  more  than  adequate  resolu¬ 
tion  in  latitude,  relative  to  the  polynomial  order  used  in  the  latitudinal 
analysis  (N=20).  In  the  cases  examined  here,  with  isolated  (in  latitude)  data 
points  and  longitude-time  data  strips,  it  makes  sense  to  extend  the  scheme  to 
latitude  as  well.  This  was  implemented  by  tying  the  maximum  latitudinal 
polynomial  order  to  the  longitudinal  wave  number  in  such  a  way  as  to  preserve 
the  modified  triangular  truncation  scheme^  throughout  the  analyses.  (Note: 

ANAINE  and  CALNE  subroutines  which  implement  this  scheme  are  included  at 
the  end  of  the  GRIDNE  program  listing  in  Appendix  C.) 
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Figure  31.  The  OOOOUT  analysis  grids  and 

jN  'j40-  gride  for  the  Stations, 
Grid  !,  and  Grid  2  analyses, 
anJ  the  Bent  OOOOUT  Ne(£40)  grid 


Two  additional  analyses  were  run  using  the  new  ,tion  scheme  -  a 

repeat  of  the  TEC  stations  analysis,  denoted  Stations  (NT),  and  an  analysis  of 
a  third  grid  data  set,  denoted  Grid  3  (NT),  which  consisted  of  data  on  an  F- 
layer  apex  grid  from  10°N  to  70°N  latitude  in  20°  increments,  0°  to  345  E 
longitude  in  15°  increments,  and  OOUT  to  24UT  in  one-hour  increments  (2880 
additional  data  points).  The  overall  statistics  from  these  analyses  are 
included  in  Table  11,  and  the  0000UT  Ng(840)  and  -Ne(340)  grids  are  shown  in 
Figure  32.  As  can  be  seen  in  the  statistics,  the  analysis  improves  slightly 
in  the  stations  case  and  dramatically  for  the  Grid  3  analysis,  which 
provides  a  better  analysis  than  the  Grid  2  case,  which  included  twice  as  much 
input  data  and  extended  down  to  the  magnetic  equator. 

A  final  series  of  limited  tests  were  carried  out  at  AFGWC  using  output 
from  the  AFGWC  TEC  preprocessor,  program  P0LY1,  and  SSIE  observations  from  the 
Astrogeophysical  Data  Base  (AGDB).  Due  to  known  deficiencies  in  the  POLYl 
program  (Bussey,  private  communication)  and  the  questionable  nature  of  the 
SSIE  data  due  to  sensor  malfunction  several  months  earlier,  only  the  effec¬ 
tive  F10  analysis  was  run  in  order  to  compare  the  mean  field  representing  each 
data  set.  The  POLYl  program  generates  24  TEC  grids,  one  per  UT ,  from  5-day 
mean  TEC  observations  and  stores  them  in  file  FOURDGRIDS.  Software  from  the 
Topside  Data  Interface  System*  was  used  to  convert  the  TEC  grids  into  grids  of 
Ne(840)  jsing  fQF2  and  M3000  from  the  ITS  ionosphere  model^  and  the  DHR 

profile  model.  (The  sunspot  number  used  in  the  ITS  model  was  36.1,  FI 0=9 1 . 2 , 

27 

which  was  obtained  from  the  effective  SSN  analysis  run  on  o-day  mean  f QF 2 
data  for  the  same  period.)  The  grids  of  Ng(840)  (51E.3  da  i  points)  were  then 
passed  to  program  GRIDNE  for  the  EF 10  analysis.  The  SSIE  N e ( 84 0 '  data  were 
taken  directly  from  the  AGDB  for  the  same  5-day  period.  This  provided  a  total 
of  4705  good  data  points  from  QMSP  satellite  F7. 

The  effective  F 1 U  values  obtained  from  the  two  analyses  were  quite 

different  -  94.71  from  the  TEC  data  and  170.61  from  the  SS I r  N  f 840 'l1  data.  As 

6 

there  was  some  question  as  to  the  absolute  accuracy  of  the  SSIE  observations 
due  to  the  sensor  malfunction,  several  additional  EHO  analysis  were  rUn  in 
which  the  SSIE  N  (840)  data  were  all  divided  by  some  constant  factor.  The 
results  are  giver  in  Table  12,  which  Msts  the  effective  F 1 0  and  average 


fp(840)  for  the  P0LY1  data  and  the  various  SSIE  tests.  It  appears  from  this 
somewhat  limited  analysis  that  the  SSIE  Ne( 840)  values  are  roughly  a  factor 
of  3.3-3. 5  too  high.  This  finding  was  discussed  with  Dr.  Fred  Rich,  AFGL/PHG, 
who  agreed  that  this  was  not  unreasonable,  as  the  potential  bias  on  the  SSIE 
sensor  on  F7  and  had  been  increased  to  its  maximum  positive  value  after  the 
March  1984  failure  of  the  ion  RPA  sensor. 


Table  12. 

Results  of  the 

SSIE  effective  F10 

analyses. 

Data  Set 

Adjustment 

Factor 

fp(840) 

EF 10 

P0LY1  TEC 

- 

1.516 

94.71 

SSIE 

1.0 

2.910 

170.61 

SSIE 

3.0 

1.680 

100.02 

SSIE 

IT 

1.642 

97.78 

SSIE 

3.3 

1.602 

95.63 

SSIE 

3.4 

1.578 

94.33 

SSIE 

3.5 

1.555 

93.08 

A  limited  amount  of  electron  sweep  data,  which  more  reliably  represent 
an  absolute  measure  of  Ng( 840)  in  the  current  sensor  status,  for  8  September 
were  compared  against  the  electron  probe  data  for  the  same  period.  A  total  of 
43  good  electro?  sweep  observations  were  found  during  the  period  examined. 
The  ratio  of  the  density  from  the  electron  probe  to  that  obtained  from  the 
electron  sweeps  ranged  from  1.5  to  5.0  with  orbit-averaged  values  of  2.7-3. 2. 
This  agrees  well  with  the  results  of  the  EF10  analysis. 


5.  40  MODEL  IMPROVEMENTS  STUDY 


The  two  tasks  which  make  up  this  study,  Tasks  1  and  3,  were  both  com¬ 
pleted  during  the  first  year  of  this  project.  The  results  were  presented  in 
detail  in  Scientific  Report  No.  1  of  this  contract^  -  only  a  short  summary 
will  be  presented  here. 

5.1  Task  1:  40  Model  Global  Analysis 

The  objective  of  this  task  was  to  review  the  AWS  40  model  program  to 
determine  what  modifications  were  required  to  configure  the  4D  model  for  a 
two-hemisphere  analysis.  The  required  changes  were  identified  and  were 
implemented  in  the  AFGWC  40  mode  code  and  documented  as  part  of  a  separate 
contract7. 

5.2  Task  3:  Model-based  4D  Height  Functions 

The  objective  of  this  task  was  to  generate  a  set  of  height  functions  for 
the  4D  model  using  a  collection  of  electron  density  profiles  generated  from 
the  I R 1 79  model.  As  it  was  decided  not  to  use  IRI79  within  the  40  model  (see 
Section  3.1),  three  sets  of  height  functions  were  generated  using  profile  sets 
constructed  from  the  DHR  model.  All  software  required  to  generate  the  model- 
based  height  profiles  and  the  associated  documentation  was  provided  to  AFGWC. 


6.  PROJECT  SUMMARY 


Over  the  two  years  of  this  contract,  a  series  of  investigations  were 
conducted  of  ways  to  expand  and  improve  the  use  of  in  situ  observations  of  the 
ionosphere  from  the  OMSP  SSIE  sensor  by  the  Air  Weather  Service  ( AWS )  40 
ionospheric  analysis  model.  The  40  analysis  model  was  modified  to  allow  it  to 
be  configured  for  a  true  global  analysis,  and  several  sets  of  40  height 
functions  were  generated  from  a  model  ionosphere  data  base.  In  the  two 
primary  investigations,  a  new  parametri2ation  of  the  topside  electron  density 
profile  and  a  technique  for  generating  northern  hemisphere  grids  of  electron 
density  from  observations  along  an  orbit  track  were  developed. 

The  topside  profile  model  developed  is  capable  of  constructing  a  topside 
electron  density  profile  using  the  full  set  of  data  available  from  the  SSIE 
sensor  (N(0+),  N(H+),  T^,  and  T  )  in  conjunction  with  ground-based  observa¬ 
tions  of  fQF2,  hmF2,  and  TEC.  An  empirical  model  of  the  height  of  the 
transition  from  0+-dominance  to  H+-dominance,  i.e.,  the  height  at  which 
N(0+)=N(H+) ,  was  developed  for  use  by  the  profile  model  should  data  from  the 
SSIE  sensor  be  unavailable. 

The  methods  developed  for  generating  electron  density  grids  from  orbit 
data,  based  on  methods  used  by  the  4D  analysis  model,  were  implemented  in  a 
computer  program,  GRIDNE,  and  tested  extensively  against  observations  gener¬ 
ated  from  a  model  ionosphere.  The  results  of  the  testing  showed  that  the 
method  reproduces  the  input  data  very  well,  but  that  the  accuracy  of  the  field 
more  than  30°  away  from  a  data  point  is  almost  solely  a  function  of  how  well 
the  Internal  model  of  the  electron  density  at  340k:n,  Ne(840),  represents  the 
shape  of  the  diurnal  and  spatial  variations  of  the  ionosphere.  The  analysis 
methods  developed  are  totally  independent  of  the  Ng(840)  model  used,  so  as 
improved  models  are  developed  they  can  be  used  to  upgrade  the  internal  Ng(84Q' 

model.  Two  internal  models  of  N  (840)  were  developed  for  use  by  program 

1  ?  ® 

GRIDNE,  one  based  on  the  Bent  ionosphere  model  and  one  based  on  the  AWS 
Damon-Hartranf t-Ramsay  (DHR)^  ionosphere  model.  The  two  models  produce 
quite  different  Ne ( 840 )  fields  for  a  given  sunspot  number,  but  without  a 
systematic  comparison  of  the  two  to  observations  it  is  impossible  to  select,  a 
priori,  the  best  of  the  two. 
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The  use  of  this  analysis  model  for  generating  joint  analyses  of  SSIE  and 
TEC  data  was  also  investigated  in  a  limited  fashion.  In  these  tests,  the  TEC 
data  were  assumed  to  have  been  converted  to  an  equivalent  Ne(840),  using  a 
profile  model,  and  then  these  Ne(840)  data  were  included  In  analyses  of 
simulated  DMSP  orbit  data.  The  tests  showed  that  this  Is  a  feasible  method 
for  combining  the  two  data  types  into  a  single  analysis. 


s 

v 

v* 


82 


R 


U1A UliLtlAZiLLlA 


REFERENCES 


1.  Secan,  O.A.,  DMSP  SSIE  Topside  Data  Interface  System,  System  Documentation, 

Report  No.  PD-NW-S2-2/5R,  Physical  Dynamics,  Inc.,  Bellevue,  WA  98009, 
September  1982. 

2.  Smiddy,  M.,  Saga  1 yn ,  R.C.,  Sullivan,  W.P.,  Wildman,  P.J.L.,  Anderson,  P.,  and 

Rich,  F.,  The  Topside  Ionosphere  Plasma  Monitor  (SSIE)  for  the  Block 
5D/F1 iqht  2  Satellite,  AFGL-TR-78-0071 ,  ADA058503,  Air  Force  Geophysics 
Laboratory,  Hanscom  AFB,  MA  01731,  March  1978. 

3.  Rich,  F.,  Smiddy,  M.,  Sagalyn,  R.C.,  Burke,  W.J.,  Anderson,  P.,  Bredesen,  S., 

and  Sullivan,  W.P.,  In-Flight  Characteristics  of  the  Topside  Ionospheric 
Monitor  (SSIE)  on  the~DMSP  Satellite  Flight  2  and  Flight  4,  AFGL-TR-80- 
0152,  ADA088879,  ATr  Force  Geophysics  Laboratory,  Hanscom  AFB,  MA 
01731,  April  1980. 

4.  Air  Force  Global  Heather  Central  Data  Format  Handbook,  Air  Force  Global 

Weather  Central ,  Offutt  AFB,  NE  681 1 3 ,  September  1983  (available 
through  Hq  Air  Weather  Service/DNTS,  Scott  AFB,  IL  62225). 

5.  Siong-Huat,  C.,  DMSP  SSIE  Uniform  Grid  Generator  (Program  SSIELD),  System 

Documentation"]  Bedford  Research” Associates,  Bedford,  MA  01730,  June 

tust: 

6.  Secan,  J.A.  and  Tascione,  T.F.,  "The  4D  Ionospheric  Objective  Analysis 

Model , "  in  Proceedings  of  the  1984  Symposium  on  *he  Effect  of  the 
Ionosphere  on  C3I  Systems']  1984. 

7.  Secan,  J.A.,  40  Ionospheric  Analysis  Model,  System  Documentation,  Report  No. 

PD-NW-83-297R-1 ,  Physical  Oynamics,  Inc.,  Bellevue,  WA  98009.  December 
1983. 

8.  Lincoln,  J.V.  and  Conkrig^t,  R.O.,  International  Reference  Ionospher~e-IRI79, 

Report  UAG-82,  World  Data  Center  for  Solar-Terrestrial  Physics,  Boulder, 
CO  80303,  November  1981. 

9.  Secan,  J.A.,  Development  of  Techniques  for  the  Use  of  DMSP  SSIE  Data  in  the 

AWS  4D  Model,  AFGL-TR-84-0167,  Aii  Force  Geophysics  Laboratory,  Hanscom 
AFB,  MA  01731,  July  1984. 

10.  Damon,  T.D.  and  Hartranft,  F.R.,  Ionospheric  Electron  Density  Profile  Model, 
Technical  Memorandum  70-3,  Aerospace  Environmental  Support  Center,  Ent 
AFB,  CO,  July  1970  (available  from  the  USAF  Air  Weather  Service  Technical 
Library,  Scott  AFB,  It  62225). 


2 

i 

j 

I  83 


11.  Flattery,  T.W.  and  Ramsay,  A.C.,  "Derivation  of  Total  Electron  Content  for 

Real  Time  Applications,"  Effect  of  the  Ionosphere  on  Space  Systems  and 
Communications.  Naval  Research  Laboratory,  Washington,  DC,  pp.  336-344, 

mv. - 

12.  Llewellyn,  S.K.  and  Bent,  R.,  Documentation  and  Description  of  BENT  Iono¬ 

spheric  Model,  AFCRL-TR-73-0657,  AD772733,  Air  ^orce  Geophsylcs  labora- 
tory,  Hanscom  AFB,  MA,  1973. 

13.  Titheridge,  J.E.,  "Ion  Transition  Heights  from  Topside  Electron  Density  Pro¬ 

files,"  Planet.  Space  Sci . ,  24,  pp.  229-245,  1976. 

14.  Bauer,  S.J.,  Physics  of  Planetary  Ionospheres,  Springer-Verlag  Publishing 

Co.,  New  York,  NY,  1973. 

15.  B’litza,  D.,  "The  Atmospheric  Explorer  C  Ionospheric  Temperatures:  Depen¬ 

dences  and  Representation,"  Report  UAS-90,  World  Data  Center  for  Solar- 
Terrestrial  Physics,  Boulder,  £6,  pp.  114-122,  May  1984. 

16.  Lockwood,  M.,  and  J.E.  Titheridge,  "Departures  from  Diffusive  Equilibrium  in 

the  Topside  F-layer  from  Satellite  Soundings,"  J.  Atmos.  Terr.  Phys.,  44, 
pp.  425-440,  1982. 

17.  Kutiev,  I.,  Heelis,  R.A.,  and  Sanatani,  S.,  "The  Behavior  of  the  0+-H*  Transi¬ 

tion  Level  at  Solar  Maximum,"  J.  Geophys.  Res.,  85,  pp.  2366-2372,  May 
1980. 

13.  Kutiev,  I.,  K.  Serafimov,  M.  Kar|dimov,  and  R. A.  Heelis,  "A  Proposed  Improve¬ 
ment  of  IRI-79  using  the  0  -H  Transition  Level,"  Report  UAG-90,  World 
Data  Center  for  Solar-Terrestrial  Physics,  Boulder,  Cu,  pp.  150-154,  May 
1984. 

19.  Titoeridge,  O.E.,  "P lasmapause  Effects  in  the  Top  Side  Ionosphere," 

Geophys.  Res.,  81,  pp.  3227-3233,  July  1976. 

20.  Moffett,  R.J.  and  S.  Quegan,  "The  Mid-Latitude  Trough  in  the  Electron  Concen¬ 

tration  of  the  Ionospheric  F-Layer:  A  Review  of  Observations  and  Model¬ 
ling,"  J.  Atmos.  Terr.  Pnys.,  45,  pp.  315-343,  May  1983. 

21.  Tascione,  T.F.,  Flattery,  T.W.,  Patterson,  V.G.,  Secan,  J.A.,  and  Taylor, 

J.W.,  "Ionospheric  Modeling  at  Air  Force  Global  Weather  Central,"  Solar- 
Terrestrial  Predictions  Proceedings,  Volume  I,  Space  Environment  Labora- 
tory,  Boulder,  CO,  pp.  367-377,  August  1979. 

22.  Thompson,  R.L.  and  Secan,  J.A.,  "Geophysical  Forecasting  at  AFGWC,"  Solar- 

-Terrestrial  Predictions  Proceedings,  Volume  I,  Space  Environment  Labor¬ 
atory,  Boulder,  CD,  pp.  350-366,  August  1979. 

23.  Automated  Data  Systems  (ADS)  Documentation,  Department  of  Defense  Standard 

7935,  February  1983. 


84 


I 


** 

/ 

I 

J 

J 

1 


24.  AFSWC  Software  Standards.  Air  Force  Global  Weather  Central,  Offutt  AFB,  NE, 

March  1983  (available  through  Hq  Air  Weather  Service/ONTS,  Scott  AFB,  ll 
62225). 

25.  Jones,  W.B.,  Graham,  R.P.,  and  leftln,  M. ,  Advances  In  Ionospheric  Map¬ 

ping  by  Numerical  Methods,  ESSA  Technical  Report  ERL  16/  -  ITS  75, 
Institude  for  Telecommunication  Sciences,  Boulder,  CO,  May  1969. 

26.  Jones,  W.B.  and  Obitts,  O.L.,  Global  Representation  of  Annual  and  Solar  Cycle 

Variation  of  foF2  Monthly  Median  1954-1956,  Ot/lTS  Research  Report  3, 
Institute  for  Telecoflinunication  Sciences,  Boulder,  CO,  October  1970. 

27.  von  Flotow,  C.S.,  "Ionospheric  Forecasting  at  Air  Force  Global  Weather  Cen¬ 

tral,"  Effect  of  the  Ionosphere  on  Space  and  Terrestrial  Systems,  Naval 
Research  Laboratory,  Washington,  DC,  pp.  385-392,  1978. 

28.  Flattery,  T.W.,  "Spectral  Models  for  Global  Analysis  and  Forecasting,"  Auto¬ 

mated  Weather  Support,  AWS-TR-242,  Air  Weather  Service,  Scott  AFBl  TT7 
pp.  42-54,  April  1971  (available  from  the  USAF  Air  Weather  Service  Tech¬ 
nical  Library,  Scott  AFB,  IL  62225). 


fi 

)i 


Y 

s 


i 

1 

1 

3 

u 

1 

3 

'i 

i 


85 


V 

A 

I 

^  .   J 


Appendix  A.  Topside  EDP  Parametrizations 


A. 1  Modified  DHR  Model 

a.  Description.  This  model  uses  a  Chapman  layer  with  a  height-varying 
scale  height  to  describe  Ng(h)  between  hmF2  and  hy,  an  exponential  layer  with 
an  exponentially  varying  scale  height  above  hy,  and  a  spline  fit  between  the 
two  layers  from  hy  -  50  km  to  hy  +  50  km. 

b.  Functional  Form. 

1 )  h  F2 <  h  <  hT  -  50 
m  i 


Ne(h)  =  Nel(h)  »  NmF2  exp  [1  -  z  -  e'2] 

h-hfflF2 
2  "  H^hT 


Hlfh)  6  Ho<VLT> 


h 

8571 


IT) 


:  empirical  DHR  scale  height  factor 


11 


:  magnetic  latitude 


IT  :  local  solar  time. 


2)  hy  -  50 <  h  <.  hy  +  50 

N  (h)  =  an  ♦  a,ch  ♦  a,.‘.h^  +  a,Ah^ 
e  vie  j 

h  =  h  -  hy 

aQ  =  R  -  25 N ' 

ay  *  0.03  (N+  -  N)  -  0.5  ( Nj  -  IT) 
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»  V-  V*  J  V-  n\ 


a2  =•  0.01  N* 


a3  =  0.0002  [(N;  -  N* )  -  0.02  (N+  -  N)] 


N  =  0.5  ( N_  +  N+) 


N'  =  0.5 


3N  3  N , 


W  *  TFT 


N.  =  Nel  (hT  -  50) 


N+  ^  Ne2  (hT  ♦  50)  • 


3)  hT  +  50<  2000 


Ne(h)  =  Ne2(h)  =  NTe‘ 


NT  -  Nel(hT) 


h^H) 


H2(h)  =  Ht  exp  (-t^) 


Ht  =  H1(hT) 


6  :  Adjustable  free  parameter. 


A. 2  Modified  Bent  Model 


l  ? 

a.  Description.  In  this  modification  of  the  Bent  model  ,  three  expo¬ 
nential  layers  are  used  to  model  Ng(h)  from  a  height  h^,  h^>hmF2,  and  a 
parabolic  layer  is  used  to  model  Ng(h)  from  hmF2  to  hj.  The  parabolic  layer 
and  the  parabolic  semi-thickness,  yt,  are  identical  to  the  Bent  model,  and  the 
scale  heights  used  in  the  exponential  layers  3re  calculated  from  the  Bent 
k-parameters. 


n 


Functional  form. 


1) 


hmF2  h  hj 


Ne(h)  *  Np(h)  -  NfflF2  I  1 


■m 


:  Bent  model  semi-thickness 


hi  s  V2  *  ^ 


Vl  *  y*  k[  -  1 


:  Bent  decay  factor  for  bottom  exponential  layer. 


2)  h,  <  h<  h-, 


Ne(h)  =  Nj  e 


h-h1 
Hi  (h) 


H^h)  =  Hq  +  a  In 


H 

a 


Hy-H 
h 
h 


In 


o 

T 

1 


H 


T 


hT-h 

N 

In  * 


I 

1 

T 


Nj  :  Density  at  hy  from  Bent  model. 


3)  hT<h<2000 


*#(h)  •  *T  • 


h-h. 


H 


2 


1 


k$  :  Bent  decay  factor  for  the  top  exponential  layer. 


A. 3  Tltherldge  M(h)  Model 

a.  Description.  Tltherldge^  developed  a  parametrized  model  for  the 
variation  of  the  plasma  scale  height  with  height  based  on  a  diffusive  equili¬ 
brium  model  of  the  topside  Ionosphere.  This  Is  used  to  define  the  variation 
of  Ne(h)  above  a  height  h^>hfflP2.  The  Bent  model  parabolic  layer  Is  used  to 
model  N^( h )  from  hfflF2  to  hj. 

b.  Functional  form. 


1) 


V2  -h<ni 
V2 


yt  :  Bent  model  semi-thickness 


hi  =  7h 


2) 


hj<  h<  2000 

h 

Ne(h)  -  Nj  exp  |  J 

ho 


=  NF2  1 
m 


*Th1s  Integration  is  carried  out  numerically  using  a  10-point  Gausian  quadra¬ 
ture  algorithm. 
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HA,hA,v  :  Mo<*e’  parameters  from  a  table*'*  based  on  h-j.  and 
the  temperature  profile  from  400  km  to  hT 
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